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Abstract 

Sea-displacement operators for fermions are defined in terms of the Fermi 
fields in a one-component Fermi system. The main conclusions of this ar- 
ticle fully corroborate the conjectures made in our earlier works, and pro- 
vide a mathematically rigorous foundation for these earlier works. These 
ideas are generalized to electron-hole systems where we are able to explore 
clearly the nature of exciton-exciton interactions. We find that exciton- 
exciton interactions in an ideal model of GaAs are not adequately treated 
simply as of the two-body type; rather the interactions are mediated by 
the exchange of other bosons that are present in this system. These bosons 
are identified explicitly and the exciton Green function is calculated. This 
exercise is also intended to be a precursor to a systematic nonperturbative 
treatment of gauge theories. 

1 Introduction 

This article is meant to address some technical issues that prevented our earlier 
work [Q from being universally embraced, and to place on a firm mathematical 
foundation some of the conjectures that appeared there. Indeed as Cune and 
Apostol H] pointed out in their very pertinent critique, there were significant 
technical drawbacks relating to the meaning of the square-root of the number 
operator in the denominator that cast doubt on the soundness of the physical 
conclusions. This article lays to rest once and for all such doubts, and as an 
added benefit, we are able to study the precise nature of exciton-exciton inter- 
actions. It is interesting to note that the rigorous formulation presented here 
is indispensable for a correct treatment of charge-conserving two-component 
Fermi systems such as excitons and relativistic electrons and positrons. This is 



1 



somewhat surprising since in our earlier work|i|| we showed that the lack of a 
proper definition of the sea-displacement operator does not invalidate the phys- 
ical conclusions in the one-component Fermi system. Here we find that a naive 
approach underestimates the nature and strength of interactions between exci- 
tons but a more careful treatment brings out many subtle features that are easy 
to overlook. Many authors who study exciton-exciton interactions assume that 
excitons are bosons that interact via two-body interactions where the exciton- 
number is conserved. Kavoulakis and Baym ||l0|| have pointed out the need to 
include Auger-like processes where the exciton-number is not conserved. Al- 
though assuming that excitons are bosons is totally acceptable(it is a matter 
of definition), there are some differences between a system such as a hydrogen 
atom (which interacts via two-body forces with other hydrogen atoms) and an 
exciton. The main difference is that in a hydrogen atom the proton and electron 
do not recombine leaving behind photons, where as an electron can recombine 
with a hole. In fact, when two excitons scatter off each other it is conceivable 
that nonradiative recombination processes take place in addition to radiative 
recombination. But in place of that, we find that an exciton can recombine 
with a special kind of electron-hole pair, which we call a solitron (to be made 
clear in the main text) and can create an intraband electron-hole pair (in other 
words, the usual kind of sea-displacement bosons found in one-component Fermi 
systems) called a valeron or conductron depending upon the nature of the band. 
Thus we find that the two-component charge-conserving undoped Fermi system 
may be thought of as consisting of several kinds of elementary excitations (which 
are postulated to be exact bosons). The usual excitons are just one class of such 
excitations. Excitons can be in a bound state or in a scattering state just like a 
hydrogen atom. Furthermore, the excitons can possess a net center-of-mass mo- 
mentum. These excitons possess an energy dispersion that is slightly different 
from the zero-center of mass excitons and are somewhat important in the anal- 
ysis. Next in order of importance we have in our system a solitron. A solitron is 
an electron-hole pair that resides at the bottom of the conduction and valence 
bands. This pair is unbound and its existence is needed in order that excitons 
interact with one another, especially when we are dealing with Auger-like pro- 
cesses where exciton number is not conserved. In addition to the solitron, we 
find the need to invoke two other kinds of bosons, valerons and conductrons. 
Valerons are intravalence band electron-hole excitations. They are analogs of 
the usual sea-bosons in one-component Fermi systems. Conductrons sim- 
ilarly are particle-hole excitations in the conduction band. All these bosons 
interact with one another and the resulting system is completely equivalent to 
the interacting Fermi system. Only the exciton couples to external radiation 
fields. Therefore, we have to consider excitons as the primary objects of interest 
and the other bosons in the system are like gauge-bosons. Material particles 
(excitons) interact by exchanging these other bosons. Furthermore, these gauge 
bosons interact amongst themselves, suggesting an analogy with nonabelian 
gauge theories. Towards the end of this article, we point out ways in which our 
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approach can be used to study gauge theories. We also compute the exciton 
Green function using the interaction terms that correspond to inelastic scatter- 
ing of the excitons off the other bosons. In future publications, we intend to 
investigate more thoroughly the practical aspects of this formalism-specifically, 
the biexciton Green function and nonlinear optical susceptibilities. 

2 Condensate Displacement Operators 

In this section, we provide some details regarding the correspondence between 
the canonical bosons and their condensate displacement analogs. These details 
are not found in our earlier works They are important since they provide 
a springboard from which we may write down the analogous statements for 
fermions. The correspondence between the sea-displacement operators and the 
canonical fermions was not made sufficiently clear in our earlier works. Indeed, 
there were inconsistencies in the technical definition, that even though had no 
impact on the physical conclusions (as we shall see in the present article once and 
for all) they did leave some room for doubting the soundness of the framework. 
Let us therefore start off with the familiar Bose systems. Let and 6k be 
canonical boson creation and annihilation operators, respectively. They obey 
the commutation rules [6k, ^j^'] — and [6k, 6j^,] — ■ Let us now introduce 
the following new objects known as condensate displacement operators, 

dq/2(q) = -^6|,6q, q^O, (1) 
V -^0 

rfo(O) = 0, (2) 

where Aq = 6960 is the number operator for the zero-momentum state. A 
word of caution regarding notation. The symbol TV" without a hat refers to a 
positive integer (a c-number) that corresponds to the total number of particles 
in the system, in both Bose as well as Fermi systems. The symbol N refers to 
the operator that corresponds to the total number of particles in the system, 
be it Bose or Fermi. The object 1 = N/N'^ is therefore the unit operator. 
The symbol Nq refers to the operator that corresponds to the total number of 
bosons in the zero-momentum state. This object does not appear in the next 
section where we deal with fermions. The square root of the operator in the 
denominator of Eq. (H) deserves special attention. In particular, when we act 
'^q/2('^) ^ state containing no particle in the zero-momentum state, we get 
an infinity multiplied by the same state, and when this further is acted on by 
60, we get a factor of zero. Zero multiplied by infinity is indeterminate. This 
tells us that the condensate displacement operator is an ill-defined operator on 
the Fock space of bosons. This is a technical problem that cannot be wished 
away. We will mitigate the severity of this problem by postulating that all 
states of the interacting system (both ground state and excited states) may 
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be expressed as linear combinations of states from a restricted Hilbert space 
that contains states of the noninteracting system with a fixed total number of 
particles, but excludes those that contain no particles in the zero-momentum 
state. Although we are unable to say precisely when this assumption breaks 
down, it is reasonable to assert that even in the case when the interactions are 
strong, either because of the intrinsic nature of the interaction or apparently 
strong due to the dimensionality of the system, the zero-momentum state of the 
interacting system will have at least one boson in it, if not a macroscopically 
large number of them. (The number operator is no longer a good quantum 
number in such systems but perhaps this is true in some average sense with small 
fluctuations around the average.) The subsequent development will justify this. 
Our assumption enables us to write down a polar decomposition of the operator 
60 as bo = exp(— iXJ)\//Vo where Xq = (1/2) (Xq + xl) is an operator that 
is strictly self-adjoint but it is not exactly but almost a canonical conjugate to 
the number operator No = b^ba > 0(in Appendix A we see that an appropriate 
interpretation of the definition of the conjugate allows us to circumvent these 
issues, for fermions we have to be more careful). On the other hand, Xq is almost 
self-adjoint but is strictly a conjugate to the number operator A'o (however, 
in the next section when we deal with fermions, Xq refers to the canonical 
conjugate of the total number of fermions). If a definition of Xq is desired, we 
may claim that it is given by the manifestly self-adjoint formula. 



(3) 



Besides the square root of the number operator in the denominator that wc have 
already made sense of, we have to make sense of the logarithm. It is defined to 
be the power-series expansion suggested by rewriting the above formula as 




XJ; = ^ In 



1+ Uo-V^o 



1 



^In 
2 




(4) 



The proof of the fact that Xq is canonically conjugate to the number operator 
Nq in the restricted Hilbert space is relegated to Appendix A. Therefore we 
may write [Xq,Nq] — i. Armed with these facts we now are able to prove 
that dq/2(q) is a canonical boson annihilation operator. We may now write 
dq/2(q) = exp(iXJ)&q. Since from the definition, exp(zXJ) depends on neither 
6q nor 6jj, we see that as far as commutation rules of dq/2(q) go, they are 
identical to those of bq since the two differ by a trivial phase that commutes 
with both these objects and their Hermitian conjugates. Next, we reproduce 
some facts that have been proved satisfactorily jlj elsewhere, and these will be 
used as the point of departure for a rigorous treatment of fermions. Earlier we 
proved [Q (for both q = and q 7^ 0), the following combined formula, 



''k+q/2^k-q/2 = A'"o^k,0'5q,0 + V A'"od-q/2 (-q)'5k+q/2,0 + %/2(q) V ^O'^k-q/2,0 
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+'^(i/2)(k+q/2)(^ + q/2)rf(i/2)(k-q/2)(k - q/2), (5) 
iVo = 7V-^4/2(q)dq/2(q). (6) 

This equation may be rewritten more elaborately and in a manner that is con- 
ducive to generalization as (q 7^ 0) 

^k+q/2^k-q/2 = V"k+q/2^k(-q) + ^k(q)V"k-q/2 

+ E <+q/2-qi/2(qi)^k-q,/2(-q + qi) (7) 

qiT^q.o 

- E ^Lq/2 + qi/2(qi)^k+qi/2(-q + qi), 



?>j,&k = ?^s(k)^ + 51 ^Lqi/2(qi)^k-qi/2(qi) 

qi#o 

- E 4+q,/2(qi)A+q,/2(qi). (8) 

qiT^o 

Here ^^(k) is the momentum distribution of noninteracting bosons at zero 
temperature, riB(k) = (5k. o^"- Also, Ak(q) = '5k-q/2,o'^q/2(q)- It is worthwhile 
to consider an alternative scheme for making sense of the definition in Eq. . 
It involves writing the number operator as 



No = iV" 



^E</2(q)'^q/2(q)| ■ (9) 



NO 

q#0 



Furthermore, we may write for the occupation number not in the zero-momentum 
state (k ^ 0) as rtk — c'(i/2)k(^)'^(i/2)k(k). With these identifications, we may 
rewrite Eq. (||) as (q 7^ 0), 

rfq/2(q) = -^li-^Y. </2(q)rfq/2(q) I bX- (10) 



qv^o 



If Eq. ( |l0|) is interpreted as a power-series expansion around unity, we can 
construct an iterative procedure to solve for (iq/2(q)- This may seem redundant 
given the fact that we have already made an elegant argument that pins down 
the meaning of dq/2(q) in terms of Xq. The reason for this new approach is 
that in the case of fermions we will not have the luxury of introducing an object 
similar to Xq for reasons that will become clear in the next section. Thus we are 
forced to seek alternatives that are more fermion-friendly. Unfortunately, these 
alternatives do not allow us to venture very far from the noninteracting case, as 
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we shall soon see. Nevertheless, this exercise is very instructive since it tells us 
that the correspondence that we write down for fermions in the next section has 
exactly the same features and are therefore correct, pending the resolution of 
the interpretation of the ubiquitous square root in the denominator. Retaining 
only the lowest order in the series expansion gives us 

rfq/2(q) = TP^O&q, (11) 

No = N"i. (12) 

Further we have (k ^ 0) ny^ = d|^^2)k(^)^(i/2)k(k) = ^{i + A^o)fek&k- But 
we also know that — 6|^6k. This seems to suggest one of two things, namely 
rik = or we should restrict ourselves to cases when Nq = >> 1. The 

latter possibility seems the most attractive until one realizes that the restricted 
Hilbert space is not so restrictive as to exclude such systems where the number 
of particles is small and finite. This will become clear in Appendix A where we 
need to only assume that the restricted Hilbert space contains no states that 
have zero particles in the zero-momentum state. Therefore we are left with the 
other possibility, rik = 0. This result is not as alarming as it seems since Eq. (^ 
is consistent only with the state where all the bosons are in the zero-momentum 
state and no bosons have higher momenta. Further iterations do not change this 
picture. That is, if we interpret the square root as a power-series expansion, we 
obtain the first order correction 

^q/2(q) « (^i + ^ E </2(qi)'^qi/2(qi) j (13) 

In order to conform to the iterative scheme, we are obliged to replace the the 
d's on the right side by the zeroth-order d's. But we know that in the zeroth 
order, d]j^^2(qi)'^qi/2(qi) = and N = N^l. Therefore the first-order d is 

the same as in the zeroth order, (iq/2(q) = ^]yo ^o^q- ^^^^ points to the 
futility of interpreting the square root as a power-series expansion around the 
noninteracting ground state. All further iterations lead to precisely the same 
result. This should tell us that one should solve the system self-consistently. 
Perhaps we should expand around an expectation value [i.e., Nq — {Nq)1 + 
{Nq — (7Vo)l)] with the expectation value determined self-consistently. This 
expectation value will be different from iV*^ even when interactions are absent, 
for example, if we consider finite temperature. Let us now try and compute 
the finite-temperature momentum distribution of noninteracting bosons. From 
elementary considerations we know that (iVo) — N^ — J2k^o(^isP^) ^'^^ ("k) = 
(^k^k)- The thermodynamic expectation values involve the chemical potential 

^"•-^ = e.p(/3(.k-M))-l- ^''^ 



6 



The reason this appears is that in the grand canonical ensemble we have to com- 
pute the trace with the Boltzmann weight (6k&k) = ^Tr{exp[—P{H—fiN)]bl^bk} 
where Z = Tr{exp[—f3{H — fiN)]} is the grand partition fimction. Using the 
cyclic property of the trace, we may write 

- ^Tr{b^e-^^"-^^^bl) 
Z 



In other words, 



Let us now try to evaluate this quantity using the condensate-displacement 



language: 



(6|,6k> = (4/2(k)dk/2(k)) 

= iTr(e-^(^-^^)4/,(k)dk/2(k)), 
H = ^ekrfJi/2)kW'^(i/2)k(k). 

k^O 



We also know that [N, <i(i/2)k(k)] = 0: 



(4,2(k)rfk/2(k)) 



lTr(e-'^(«-^^)4/2(k)4/2(k)) 



Z 

= lTr(e-''(^-^^)e''(^-^^'dk/2(k)e-^(^-''^)4/2(k)), 

(4/2(k)rfk/2(k)) = e-^^'^lTr(e-'^(«-^^)4/2(k)4/2(k)) 
= e-^''^(l + (4,2(k)4/2(k))). 



(4/2(k)dk/2(k)> - — nr^^- 



The crucial /i seems to be missing. The reason is somewhat subtle [Q. It has to 
do with the fact that the trace in the original case using the parent bosons spans 
all states including those with Nq — 0. However, in the condensate-displacement 
language, the trace is over all states except those that have Nq — 0. There are 
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quite a number of states that have A^o = and an arbitrary total number of 
particles. It is therefore not surprising that we have encountered a discrepancy. 
The best way to resolve this is to introduce a Lagrange multiplier that allows 
us to control how many bosons there are in states with zero momentum. Thus 
when it comes to taking the trace over states in the condensate-displacement 
language, we have to be careful to include a new chemical potential that couples 
to Nq and not just to N. When this is done we can easily show 

These considerations carry over to fermions where the difhculties are under- 
standably far more severe. After all, trying to describe fermions using Bose-like 
objects, and to do it exactly, is a daunting task. We shall now use this insight 
to write down a correspondence for fermions. 



3 The Nature of Sea-Displacement Operators 

In this section, we provide details of the correspondence between the sea-displacement 
operators and canonical fermions. The sea-displacement operators have been in- 
troduced elsewhereQ but an explicit formula for the sea-displacement operator 
in terms of the Fermi fields was lacking. In this section, we once and for all 
pin down the definition of the sea-displacement operator and show that they 
obey commutation rules that are somewhat complicated but in the limit of a 
generalized random-phase approximation (RPA) (to be made clear in Appendix 
B) they obey canonical boson commutation rules. This exercise hopefully is the 
final word as far as the technicalities of the correspondence goes. Furthermore, 
as we shall see in the conclusion, the rigorous formulation reinforces the rather 
startling claim[Q that in ID, so long as the interaction between the fermions is 
purely repulsive and possesses a Fourier transform, and is sufficiently weak so 
that all states of the interacting system (both ground state as well as excited 
states) may be expressed as linear combinations of low lying excited states of 
the noninteracting system and its ground state, then the momentum distribu- 
tion of the interacting system possesses a sharp Fermi surface at precisely the 
same place as the noninteracting system. Let us now proceed to the main task 
at hand. 

Let Ck and be canonical fermion annihilation and creation operators. We 
may write {ck, c^^'} — and {ck, c|^, } = f^k k' • "^k = ^j^^k denote the number 
operator. The sea-displacement operators for fermions are defined as (q ^ 0) 

^k(q) = np(k-q/2)[l-nf(k + q/2)] ^ 4-q/2gk+q/2 (15) 

Ak(0) = (16) 
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where ufO^) — ^(fc_F— |k|) is the zero-temperature Fermi distribution. As before, 
the square root of the number operator in the denominator in the definition 
requires clarification, a point noted by Cune and Apostol[||. Their approach 
for deahng with this problem is unfortunately, not adequate!^. Neither it 
seems is the approach used for dealing with this problem in the case of bosons, 
namely that we be able to interpret the object U(k) = -^^c^^ as being a unitary 
operator. This reason for this additional complication probably stems from the 
fact that in order for this object to be unitary, we should be able to find a self- 
adjoint canonical conjugate Pk of the number operator such that [-Pk,'^k'] ~ 
iS^ k' • Such an object is not likely to exist not only because of the positivity 
of the number operator nk, but also because of idempotence [Q. In fact, 
it is a simple matter to convince ourselves that idempotence is inconsistent 
with the existence of Pk. On the one hand, we have [Pk,'T-k] = [Pk, 'T-kl'T-k + 
'^klPkj'^k] = Sink. On the other hand, [Pk,n^] = [Pk,'^k] = *• This suggests 
that nk = (l/2)l-a state of affairs seldom realized if at all. Some readers of this 
work familiar with the more traditional bosonization approaches may point to 
the importance of 'point-splitting'-'a procedure that does not allow us to write 
TZk = cj^Ck. They may wish to suggest this as the main reason for all these 
difficulties. It is likely that point-splitting is a necessary technical consideration 
only in systems that have kp oo and m —^ oo but vp = kp/m < oo, in 
other words, a system where both the mass of the particle and the density are 
infinite such that the Fermi velocity is finite We deal with systems that 

are more physical (finite kp and m) and are, by and large, immune to these 
considerations. While the metaphorical definition of Eq.( |l^) is quite adequate 
for most practical purposes, it is desirable to have a better understanding of the 
meaning of the square root. 

Let us divert our attention for some time to the more mechanical aspects 
of this program in the hope that soon we will be able to address this delicate 
technical problem of the square root in the denominator. In a manner entirely 
analogous to that of bosons we may write 

4+q/2Ck-q/2 = V"k+q/2-4k(-q) + Ak(q) V"k-q/2 

+^ E ^Uq/2-qi/2(qi)A-qi/2(-q + qi) 

qi7^q,o 

E ^Lq/2 + qi/2(qi)A+q,/2(-q+qi). (18) 

qiT^q.o 

Furthermore, 

nk = = np{k)^+\ ^Lqi/2(qi)^k-qi/2(qi)-A J2 ^Ic+qi/2(qi)^k+qi/2(qi) 

qi#0 qi#0 

(19) 

where np{\i) — 9{kp — |k|) is the zero-temperature Fermi distribution function. 
The parameter A is a book-keeping device, familiar from perturbation theory. 
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where in the end we set A = 1. This is a very useful device as we shall soon 
find out, one that partially addresses the problem of the square root in the 
denominator. The proof of Eqs. (|l^) and ( ]l9| ) rests crucially on the meaning of 
the definition of Ai^{q) shown in Eq. (p5|). Let us now consider the problem of 
ascertaining the meaning of Ak(q) in a 'perturbative' manner. That is, we shall 
interpret the square root of the number operator as being the square root of the 
Eq. (|l|) where we find a natural unit operator about which we may expand the 
square root. Thus if we consider the lowest order in A we have rik = np{'k)l. 
Since ni(;-q/2 = 1 when np(k — q/2) = 1, we may write, 

A(q) = "F(k - q/2)[l - npik + <l/2)]cl_^^^c^+^/2. (20) 



These equations obviously solve Eqs. ( p^ ) and (19) in the lowest order A" = 
1 except that we seem to have the additional curiosity c|^_|.q/2''k-q/2 = if 
n_F(k+q/2) = 1 and n_F(k-q/2) = 1 or if n_F(k+q/2) = and n_F(k-q/2) = 0. 
This is hardly surprising given the fact that — np{\i)l is consistent with our 
restricted Hilbert space having exactly one element, namely the ground state 
of the noninteracting Fermi sea ji^S*). This means cjj.^q^2'^k-q/2l-^'S') = when 
n_F(k+q/2) = 1 and n_F(k-q/2) = 1 or if n_F(k+q/2) = and n_F(k-q/2) = 0, 
a fact easily verified. Parenthetically, we note that Eq. was anticipated in 
the comment by Cune and Apostol Just as in the case of bosons, further 
iterations do not affect these conclusions. Therefore, if one starts with the 
noninteracting system and expands around the noninteracting ground state, one 
is able to treat only the noninteracting system. Just for the sake of completeness, 
let us see where the first order terms lead us. When rii?(k+q/2) — 1 and rii?(k— 
q/2) = 0, we have from the definition Eq. ( p^ ) a result that is independent of 
A and therefore valid to all orders in A, Y^7ik+q72^k(— q) = '^t.+q/2'^^-i/^- 
addition, we have ^k(q)\/'T-k-q/2 ~ cj^-|_q/9Ck-q/2 when np(k + q/2) = and 
^^(k — q/2) = 1. These two identifications are entirely consistent with Eq. 
( p^l ) to all orders in A. The other cases namely, ni?(k ± q/2) = and npCkzt 
q/2) = 1 are given below. They have to be proven iteratively. When npCk + 
q/2) = and np{k — q/2) = 0, we can see from Eq. (|7|) that cjj.^q^2'^k-q/2 = 

-^Eqi5^q,o^k+q/2-qi/2(qi)^k-qi/2(-q + qi)- Wc uow show that this rcsult is 
identically zero. This is consistent with our earlier claim that further iterations 
do not change the zeroth order results. Since, now if we replace the A's by the 
zeroth order ones, we get a first order result for c|j.^qy2'^k-q/2 equal to zero: 

clc+q/2^k-q/2 = A ^ 4+q/2-qi/2(qi)^k^qi/2(-q + qi) 

= A ^ '^f(k + q/2-qi)4^q/2Ck+q/2-qi4+q/2_q,Ck-q/2 
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= A 



npik + q/2 - qi)(i - nk+q/2-qj4+q/2Ck-q/2- 



qi^q.O 

The number operator ?^k+q/2-qi — «F(k + q/2 — qi)l, and therefore the first 
order result is zero as welL Similarly, we find the same answer when ^^(k + 
q/2) = 1 and np^k — q/2) = 1. To put it another way, iterations around 
nk = rii? (k)l do not change the form of nk- That is, nk remains frozen at the 
value riF (k)l. This means that the only state consistent with such an identity 
is the noninteracting ground state. Using this, we can convince ourselves that 
the program is consistent, just as the program was consistent in the Bose case 
which has been proved rigorously by other means. 

In order to do better than just remain at the noninteracting ground state, 
one must, just as in the Bose case replace nk = (^k)! + (fi-k — (nk)l), and 
expand around an expectation value (rik) 7^ for all k. Thus we are obliged 
to consider a system with interactions and at a finite temperature in order for 
this scheme to be of practical significance. The philosophy is that we solve 
for (rik) self-consistently and then pass to the limit of weak interactions or low 
temperature if one wants to study the ideal case. Therefore, we may write 



1 1 



V'^k-q/2 \/("k-q/2) 



^ nk-q/2 - (?^k-q/2)^ 
("k-q/2) 



(21) 



Let us now use the insight obtained in the case of bosons to write down the 
correct commutation rules obeyed by the sea-displacement operators. This is a 
nontrivial task given the subtleties involved. The Bose case was soluble exactly 
via a polar decomposition which we do not have here. We write down here 
the final answers, details of which may be found in Appendix B. The exact 
commutation rules obeyed by the sea-displacement operators it seems are rather 
complicated. The sea-displacement operators do, however, obey exact closed 
commutation rules as we shall see in Appendix B. There is a natural and simple 
regime where the approximate commutation rules are those of canonical bosons. 
It is the regime of the RPA of Bohm and Pines Q . They are 

Ak(q),Ak'(q')l_^ = 0, (22) 

J RPA 

Ak(q),4,(q')]^^^ = n^^(k~q/2)[l-n;^(k + q/2)]<5k,k"5q,q'. (23) 

We saw in the Bose case that the RPA-commutation rules may be taken to be the 
exact rules the only care that we had to exercise was to couple the condensate 
bosons to their own chemical potential. Here too we have to couple the sea- 
bosons to a momentum dependent chemical potential as we will soon show. We 
find that it is quite appropriate to view the objects Ak(q) as being exact bosons, 
although for special values of the indices they behave rather strangely. These 
simple-looking results belie the notoriously difficult and technical problem of 
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the square root in the denominator, which incidentally, we have only partially 
resolved [see Eq. (^)]. One potential criticism of this is to claim that a scheme 
such as Eq. is inconsistent with the exact commutation rules presented in 
Appendix B. The rebuttal to such a critique is that even Eq. ( |2l| ) is approximate, 
no matter how many orders are summed-the reason being that this assumes that 
fluctuations in the number operator are small compared with the mean. This 
is inconsistent since we may show for example that (n^) — (rik)^ = {n-k}{l — 
(rik)). Therefore, the more noideal the momentum distribution is, the more it 
fluctuates. In fact, a momentum distribution that is highly nonideal (nt) ~ 1/2 
for most momenta, has the largest fluctuation equal to the mean itself. There 
is another important point that should be mentioned. The definition in Eq. 
(^5|) closely resembles the definition of condensate displacement operators for 
bosons. This is certainly a plus. Further, the factor ni?(k— q/2)[l — nF(k+q/2)] 
serves as a cutoff function. It fulfills a very important role. It ensures that 
the kinetic energy operator is positive-definite. The authors initially tried a 
number of alternatives that attempted to include terms beyond the RPA-like 
term, for instance one with [1 — ni?(k + q/2)] [1 — ni7(k — q/2)]. Such attempts 
are always unsuccessful for the simple reason that they lead to non-positive 
kinetic energy operators. With the cutoff function we have introduced, we may 
write the kinetic energy operator in the sea-displacement language as K = 
fco^ + I]k,q7^o^k(q)-4j^(q)^k(q)- it can be seen that Wk(q) = (k • q/m)nF(k- 
q/2)[l — np{h + q/2)] > 0. Here fco = X^k £k'^_F(k) is the kinetic energy per 
particle. In order not to disappoint the attentive reader, we collect here some 
facts that are true in the absolute sense. From the exact definition in Eq. ( p^ 
it is fairly obvious that (independent of the meaning of the square root) 

[^k(q), np] = Ak(q)(5p^k+q/2 - '5p,k-q/2)- (24) 

More generally, 

Fi[{np + <5p,k+q/2 - '^P,k-q/2}])Ak(q) = A^{q)F{[n]) (25) 
From the definition in Eq.( p^ ) we may write, 

Ak(q)A[(q) = np{k - q/2)[l - npik + q/2)](l - nk+q/2) (26) 

The other commutation rules are collected in Appendix B. They are important 
only for special values of the indices k, k , q and q . In the RPA-sense and 
whenever we are willing to overlook these nuances (which is always, practically 
speaking) the RPA commutation rules of Eq. ( |2|) and Eq. ( |2^) are quite suffi- 
cient and easy to use. Let us now try and compute the properties of the free 
Fermi theory using the machinery that we have just laid down. For this we need 
to first ask : How to express A|^(q)Ak(q)? We expect the answer to depend 
crucially on the meaning of the square root. We do not have any more in- 
sight into the nature of this object, but we will be needing its expectation value 
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when we try to compute the finite temperature properties of noninteracting 
Fermi system in the sea-displacement language. Fortunately, the expectation 
value is obtained by a reasonably straightforward method and one that is very 
reminiscent of the Bose case. We shall have occasion to introduce a chemical 
potential that is momentum dependent and scales as the logarithm of the vol- 
ume of the system so that one may venture into regions where the square root in 
the denominator vanishes (as does the numerator) where the sea-displacement 
method breaks down. The introduction of the chemical potential restores the 
meaning of the trace that in the Fermi language was intended to span over 
all the states. Let us now evaluate the finite-temperature properties, specifi- 
cally, the finite-temperature momentum distribution of noninteracting fermions 
in the sea-displacement language. To this end, let us now compute the following 
correlation function: 

Z(k-q/2) = Tr(e-'^(^-^^-^''-i/^"''-<i/2)). (28) 
Using the cyclic property of the trace, we have 

(4(q)A(q)) = z{\i-ci/2f '' [A(q)e-''^^-'''^-^''-/="'^-''/=)4(q) 



' -Tr 



-/3(-H-A'W-Mk-q/2nk-q/2)p/3(-ff-M-'V-Mk-q/2nk-q/2) 



(29) 



Z(k - q/2) 

e'/3i^g-/3Mk-q/2 ^ q/2) ^'' [e"'^*'^"^^"^''"''^""'"''^'^^k(q)4(q) 



^-<./^(Ak(q)4(q)). 
But we know from Eq. ( p6| ) that 

(Ak(q)4(q)) = n;^(k-q/2)[l-nj.(k + q/2)][l-n'5(k + q/2)] 

« riF(k-q/2)[l-nj.(k + q/2)] (30) 

at low temperatures, since the thermodynamic expectation value n''(k -I- q/2) 
vanishes exponentially when ni?(k-|-q/2) = 0. Therefore we may write (A|^(q)Ak(q)) 
ni?(k — q/2)[l — np{k -|- q/2)] exp(— /3k • q/m) exp(— /3^k-q/2)- The chemical 
potential /ik has to be chosen so that the correct thermodynamic expectation 
values are recovered. It is clear that in the absence of the chemical potential, we 
are bound to obtain answers that are incorrect, unless one is at absolute zero. It 
is also clear that in order for Eq. (^l|) to make sense, one must choose a chemical 
potential that scales as the logarithm of the volume of the system. Let us set 
exp(/3^k) = l^exp(Ak) for some Ak. This may be determined by demanding the 
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self-consistency of 



q#0 

? kg 



J2 ^F(k)(l - 7iF(k + q))e-^^e-'^'<'e-'^''^ (31) 



1^0 

1 —e^^ Z2^^^ nFiq)]e- 

qT^k 



+ [l-nf(fc)]e-^^i^n^.((Z)e'5^^e-^«. (32) 



V 
q 



Let us now focus on a system in two space dimensions, in which case 
nFik)e-^' = f^\nF{k){l~nFAk))e-^^ef'& 



/2^\ 
\ m J 



nF{k), 

L^„Me-i^^, . ^£'2,,d,(|^)(l-„„(,)).^S 

The last result follows if we assume, as we shall, that we are working at low 
enough temperatures ksT « Ef- For this program to be self-consistent, we 
have to make this assumption. Perhaps there is a way out of it, but for now the 
authors are unable to find it. In any event, it is this regime that is of considerable 
physical significance and the fact that the sea-displacement scheme is consistent 
in this regime is very comforting indeed. Moving on, we write 

lEll-M*-*^^(^)e-&^ (33, 

Substituting these back into Eq. ( |3^ ) we find nF^pik) = riF(k)ni?_^(k) -|- [1 — 
ni^(k)] exp[— /3(fc^ — ky)/2m]. At low temperatures when |k| > kp^ we know 
that UF^fiik) ~ exp[— /3(A;^ — fc|.)/2m]. Therefore, the sea-displacement lan- 
guage gives the exact same result as the original Fermi language. Having suc- 
cessfully fixed the chemical potential of the sea-displacements, we now turn 
to the problem of computing the correlation functions of the noninteracting 
system at finite temperature. For instance we have the four-point function, 
-F(kq; k'q') = (Ck-Hq/2Ck-q/2cj^'_q'/2'^k'+q' /2)- ^'^^ Fermi language, it is sim- 
ply given by 

F(kq;kq') -5k,k"5q,q'«F,Mk + q/2)[l-nj.,^(k-q/2)] (34) 
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It may be shown that in the sea-displacement language, the same quantity may 
be evaluated as 

F(kq;k'q') = 

^k,k"5q,q' [(V"k+q/2^k(-q)^j^-(-q')^nk'+q72) + V"k-q/2 V'^k' -q72^k' 

+ '^k,k' ^q,q' (4+q/2-qi/2 (^1 ) A+q/2-qi /2 (QI )) (k+q/2-qi ) [l-TZf (k-q/2)] 

qiT^q.o 

+ '5k,k"5q,q' J2 K-q/2+qi/2(qi)A-q/2+qi/2(qi))'^i=^(k+q/2)[l-«i=^(k-q/2 + qi)]. 

qiT^q^o 

_ _ (35) 

Equation ( |35| ) is obtained by placing two expressions in Eq. (|l^) next to each 
other, contracting the indices and finally using Eq. (^^. Insofar as the scheme 
for computing the momentum distribution is satisfactory, it may shown that 
formula in Eq. ( ^ ) is in fact identical to (^). To show this rigorously, we have to 

properly interpret the average (^nk+q/2^k(— q)^|^' (— q )\/"-k'+q72)- ^^'^^ 
is a simple matter to evaluate this quantity exactly in the Fermi language. From 
Eq. (H) we have ^nk+q/2-4k(-q) = c^q/2'^k-q/2'^i=^(k + q/2)[l-ni;^ (k-q/2)]. 
Therefore we may write quite unambiguously 

(V'^k+q/2^k(-q)^J^. (-q ) ^«k'+q72> 
= ufA^ + q/2)[l - npAk - <l/2)]npik + q/2)[l - npik - (l/2)]5^ySl^^) 

Similarly we may write, 

(^lc(q)V'^k-q/2^?v-q72^k' (q')> 

= 72fA^ + q/2)[l - ufA^ - <l/2)]nFik - q/2)[l - nF{k + q/2)],5k,k"5^3?) 

The main conclusion of these arguments is that one must be careful when eval- 
uating these expectation values. In particular, it is wrong to make the approx- 
imation Tik+q/2 ~ (^k+q/2)lj slucc this implics 

(^Uq) V^hl^^n^' -q72 A' (q ) > =? nf^(k-q/2)nj.(k'-q72)(4(q)Ak'(q')) 

= n^(k-q/2)[l-n^(k + q/2)]Vk"5q,q' ^'^"^ (^) ■ (38) 

This result is inconsistent with the correct result in Eq. that does not vanish 
in the thermodynamic limit. Having said this, it is still not a poor approxima- 
tion at low temperatures (fcsT << Ep)- The reason being that although the 
approximation on the right-hand side of Eq. ( ^8| ) vanishes in the thermodynamic 
limit at any finite temperature, the exact answer for this correlation function 
as derived above also vanishes for a different reason. At low temperatures, 
due to the presence of a product such as np^pik + q/2)[l — ripik + q/2)], the 
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correlation function vanishes exponentially fast as exp[—{ey^^q/2 — EF)/kBT]. 
Therefore, while it is certainly advisable to be cautious, we are allowed some 
leeway at low temperatures. Having done all this, it is now a simple matter to 
convince ourselves that Eq. ( |3^ ) is identical to ( p^). If ni?(k + q/2) = 1 and 
ni;'(k— q/2) = 0, then all but the first term in Eq. ( |35|) vanish and the identity is 
trivially satisfied. If rii?(k + q/2) — and np(k~q/2) — 1, then all but the sec- 
ond term vanish and the identity is again trivially satisfied. What remains now 
is to compute the terms when ni;'(k + q/2) = and ni?(k — q/2) = and when 
nF(k-|-q/2) = 1 and ni?(k— q/2) = 1. In the first instance, we have to evaluate 
the sum = Eqi^q,o(^k+q/2-qi/2(qi)^k+q/2-qi/2(qi))nF(k+q/2-qi). For 

this we have to make use of the simplified expression in 2D (Aj^(q)Ai(;(q)) = 
ni?(k — q/2)(l — npi\L + q/2)) cxp(— /?k • q/m)27r/3/ {mV). Therefore, we may 
write 5*1 as 

At low temperatures we may approximate Eq. ( ^4| ) as -F(kq; kq) = exp{/3[fc|, — 
(k+q/2)2]/2m} since when nF(k-q/2) = and Ep » kT, and n_F(k+q/2) = 
implies |k + q/2| > kp. Therefore we have 1 — np^pCk — q/2) w 1 and 
np.pCk + q/2) — exp{(3[kp — (k + q/2)^]/2m}, and hence the result follows. 
Similarly, when ni?(k + q/2) — 1 and npCk — q/2) = 1 we have 

H (^lc-q/2+qi/2(qi)^k-q/2+qi/2(qi)>[l " np{k - q/2 + qi)] 
qi^q.o 

= 77 2^(1 -f^F(Ql))e ( — 1 e'^ = . (40) 

Qi V m y 

Analogous to the earlier case, we may approximate Eq. (|3^ ) as -F'(kq; kq) = 
exp{/3[(k — q/2)^ — kp]/2m}. Therefore, this scheme gives the right answers 
at finite temperature for both the momentum distribution and the four-point 
functions. 



3.1 Role of Partial Isometries 

The term partial isometry |^ refers to, intuitively, an almost unitary operator. 
The absolute unitary nature is spoiled by an object that is by now ubiquitous 
in this article, namely the canonical conjugate Xq of the total number operator. 
Thus we would like to interpret the object 

U{k)^np{k)^cl (41) 

as an almost unitary operator when rti?(k) ~ 1. The presence of the square root 
of the number operator in the denominator led others O] to conclude that our 
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theory contains divergences. A closer examination tells us that this is not the 
case. The creation operator cj^ scales as the square root of the number operator. 
In other words, if we write 

4 = V^exp{-ie([n];k)}exp(iP^) (42) 

, the creation operator becomes a nontrivial complex operator of unit modulus 
(ideally) times the square root of the number operator. The creation operator 
differs from the square root of the number operator by a phase factor. Analogous 
to what we found in our earlier works ^ , this phase factor consists of a canonical 
conjugate to the number operator denoted by Pk and a real functional of the 
number operator 0([n];k). The sea-displacement operator may be written in 
the following suggestive form: 

Ak(q) - np{k - q/2)[l - npik + q/2)]e-'«(["l^''-i/2)e' (43) 

This form, although more manageable than the definition in Eq. (p^), is am- 
biguous since we have yet to pin down the meaning of Pk and 0([n];k). We 
shall adopt an approach similar to the one we outlined in our earlier works 
where we transform to the Fourier space and write down a formula for Pk in 
terms of so-called momentum currents (as opposed to real currents). As before, 
we find in our formalism the object Xq that spoils the unitary nature of these 
operators. This object though troublesome, is very important since without 
it our program would be inconsistent. To understand the need for this object 
more clearly let us write (as we rightly anticipate) Pk = Xq + Pk where Pk 
is manifestly self-adjoint whereas Xq is not self-adjoint and obeys the relation 
[Xq, N] — i. If Xq is interpreted as being self adjoint, then the resulting formula 
for A|^(q)Ak(q) obtained from Eq. ( ^ ) would be inconsistent with Eq. ( p^ 
unless nk = nFO<-)i, a state of affairs realized only in the noninteracting system 
at zero temperature. Therefore, it is imperative that we do not interpret Xq 
as being self-adjoint. However, just as in the Bose case we may still choose the 
hermitian part of Xq provided we interpret the polar decomposition in Eq. ( ^2|) 
in the distribution theoretic sense to be made clear below. Furthermore, it is 
not clear what the role of idempotence is in this formalism. In Appendix B, we 
write down closed commutation rules obeyed by the sea-displacement operators. 
These are the exact rules. The derivation of these do not require the use of idem- 
potence, but we have to make use of the fact (ck)^ = 0. This and other similar 
observations such as idempotence require point-splitting techniques. This pro- 
cedure is well-known to the traditional bosonizing community, mainly because 
it is indispensable in dealing with relativistic systems in 1-f-l dimensions which 
have a different mathematical structure than the nonrelativistic systems that 
we consider here. Perhaps these techniques can even be made mathematically 
rigorous in those systems. We had to consider the field operators in real space 
to be operator- valued distributions [Q, which required the introduction of a sep- 
arable one-particle Hilbert space, which amounts to placing the system in a box 
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with periodic boundary conditions. Similarly, we are now forced, it seems, to 
do the same with momentum space, that is, we postulate in addition a short 
distance cut-off that amounts to introducing a lattice constant. Space is now 
discretized and finite as well. There is a large macroscopic length scale L (di- 
mension of the box) and a small microscopic length scale a (the lattice spacing) . 
In the end we set L oo and a — *■ and hope that our edifice remains intact. 
Let us try and make more sense out of Eq.( |43|) . First, we would like a formula 
for Xq in terms of the Fermi fields without involving 0([n];k). To this end, let 
us postulate the existence of a 'super boson'. It is defined as, 

-iX, ' 



C^e-'^^'WN (44) 

where N = X^k'^k'^k and [Xo,iV] = i. If Xq is strictly self-adjoint, then C is 
indeed an exact boson. Otherwise it is a quasi-boson. We would not like to have 
Xq as self-adjoint for reasons already alluded to. The claim is that the filled 
Fermi sea is obtained by creating TV" number of these bosons from the vacuum 
of the fermions. 

iC^r"\0) = \FS)= I n 41 |0) (45) 

\|k|<fei=. 

We may rewrite the filled Fermi sea as follows. 



1^^) = ( n 4) |0) =exp(^n^^(k)?n(4) ) |0) 

VlkKfej. / \ k / 



(46) 



Jk|<fej. 

From this we may deduce. 



= exp (^-^ "F(k)Zn (4) ^ w ^e'^o (47) 

upto some phases that we will soon show are unimportant. From this we may 
read off a formula for Xq. 

^" = ]^E"^W'"(^k) (48) 

k 

From this it is clear that Xq is not self-adjoint and further we have to show that 
it is a conjugate to N. For this let us compute, 

c'^^Xqc-'^'^ =Xq + X (49) 
Using the definition in Eq.( Q) we have, 

e^^^Xoe-^'^ = ^Y1 "^C')'" (e'^^Cke-^^*) = Xq + X (50) 
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Therefore the conjugate obeys aU the expected properties, independent of the 
meaning of the logarithm. We would now like to ascertain how this object 
commutes with other objects such as Up = cJ,Cp. 

e^^"pXoe-'^"p = ^Y1 ^^-C^)'" (e^'^-Cke-'^"") = + (51) 

k 

Now we would like a formula for the self-adjoint operator Pk- First we would 
like to point out some obvious pitfalls. If one takes the point of view advocated 
in our earlier work||l[, then we have to face up to the fact that the number 
operator obeys idcmpotence since only the real space quantities are distributions 
whereas the momentum space quantities are bona fide operators. Idcmpotence, 
unfortunately is inconsistent with the existence of a canonical conjugate such 
as Pk as has been argued earlier. This then means that somehow we have to 
eschew idempotence in favor of retaining the conjugate. It is still unclear to 
the authors why idempotence is such a big hurdle and how to overcome it. We 
take the point of view that any Pk should result in a Ak(q) that obeys Eq.( pj[) 
and furthermore we also demand as indicated earlier, the identity in Eq.( |26| ). 
This matter has been studied by us quite thoroughly and we are unable to find 
a completely satisfactory answer. However, we write down one possibility. This 
definition involves the problematic line integral encountered in our earlier work 
as well[|l|. The resolution to this may be obtained by transforming to the real 
space analogous to the approach used for making sense of the DPVAll). The 
quantity Pk may be written down as follows. 

Pk=y" dl{~l/np)I{p)+e{[n]-k)- J dl [-ie,VP](p) (52) 

The line integral is to be interpreted as being carried out after transforming 
the quantities to real space analogous to the DPVA||l|. The quantity I(p) = 
(l/2z)[c|,Vcp — (Vc|,)cp]. Finally, we have to write down a prescription for Q. 
Since we have discretised the momenta by placing the system in a box, we may 
use the following natural ansatz. There is no need to make contact with the 
free theory very likely since in the number-phase respresentation the free theory 
is recovered automatically. Let kn — (27r/P)(rii, n2, na). Define the mapping 
t(kn) = 2"i X 3"^ X 5"^. It may be seen that i is a bijection. From this we may 
write [p^, 

e(H;k„) = 7r^0(t(k„)-t(p™))np„ (53) 

Pn, 

Here 9{x) = 1 if a; > and 6{x) = is a; < is the usual Heaviside step function. 
This obeys the required recursion relation below, necessary for ensuring that 
Fermi statistics are obeyed. 

0(["P1 - ^p,,k'}];k) + e([n];k') - e(N; k) - e([{np, - 5p,k}];k') = ±7T (54) 
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3.2 Short-Range Interactions 

We have found a curious feature in our investigations of the nature of the RPA 
that is worth mentioning. It seems that depending upon how one groups the 
Fermi field operators, one encounters an enigmatic duahty between repulsion 
and attraction(see Ref. 0, chapter 4). It is as follows. Consider the interaction 
termjl^ in jeUium j^: 

^' = E 1^ E 4+q/2^k-q/2cj^._q/2Ck'+q/2, (55) 

\ " t t 

2V ^ '^k+q/2Ck'+q/2V-q/2'^k-q/2 

= " 51 XI ^2v' "^P+Q/2^P-Q/2^P'-Q/2'^P'+Q/2- (^^) 
p#p' Q#o 

The last result is obtained simply by relabeling the indices. If we set (see end of 
Appendix B) c^qy2Ck-q/2 =rpa [Ak(-q) + ^k(q)] with equality to the level 
of RPA, then we may write 

= E 1^ E [^k(-q) + 4(q)]K' (q) + 4 (-q)]- (57) 

q^^o k^k' 
Alternatively, we may also write 

^/ = - E E ^[^k(-q) + 4(q)][A' (q) + 4 (-q)]- (^8) 

The first alternative yields results identical to the ones presented in our earlier 
work Q . The momentum distributions derived from this interaction term pos- 
sess sharp discontinuities. The second alternative is more interesting. It has 
a minus sign signifying an apparent change from repulsion to attraction (only 
apparently so). It has been argued in the literature[0 that if we assume that 
^^q = ^'^o is independent of q (in other words, a (5- function interaction in real 
space), then we may conclude by examining Eqs. ( |55|) and (^6|), that Hj = —Hj. 
Therefore, Hi = 0. The exchange correlation energy of a system of fermions 
interacting via a S -function interaction is identically zero. This is not surprising 
since we know that two fermions cannot be at the same point and the interac- 
tion is zero unless they are at the same point, therefore one obtains an answer 
zero. This means that for such systems the Hamiltonian is a simple functional 
of the number operator. 

H = fk^k - 1^ E "k+q/2'^k-q/2- (59) 

k k,q5^0 
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Unfortunately, this Hamiltonian does not yield results too different from the 
noninteracting case. However, we have found that if one relaxes the assumption 
that Vq is strictly independent of q and instead say, puts in a weak dependence 
on q, it is still a good approximation to ignore the exchange correlation energy 
(at least in comparison with the kinetic energy and the exchange self-energy), 
and we obtain some rather interesting results for the one-particle Green function. 
For instance, we could assume that the potential in real space instead of being 
a strict (5-function is somewhat smeared. For example in ID, we might take 
V{x) — e'^ {a / tt) {x'^ + a^)~^. We merely quote the final results for the one 
particle properties as this is just an aside. In order to derive these results one 
has to make use of the Schwinger's functional approach to dealing with quantum 
field theories |]. This exercise highlights the importance of fluctuations in the 
momentum distribution in determining the salient features of the one-particle 
Green function in systems that have very short range interactions and are not 
too strong. Following the notation of Kadanoff and Baym we write down 
the spectral function and the spectral width as 



AM . (eo) 

r(p,c^) = ^/-k{p,lu) (61) 

where k(p,uj) ~ {ui — ip + /i)^ — 4i^(p) < 0. In the event k{p,uj) > 0, 
\/—k{p, oj) — 0, by definition. Here i^(k) is related to the number-number 
correlation function A^(k, k ) = {n\^n-^^') — {n]^{n-^^i) where ek is the modified 
dispersion relation. Ck — fk — ^^^^ X^q^^o "q('*k-q)- Furthermore, F(k) = 

qVo^qV^('* - q, k - q ). The momentum distribution has to be 
determined self-consistently. 

(v^fi)rd« — — ^ (62) 

J J -7^/2 ^f3(i^- t.)f,2(Sy/F(^) sine 

These formulas are incomplete unless one has a prescription for the object 
A^(k, k). We note here some properties of this object. First, iV(k, k) = 
(nk)(l — {n\s))- It may occur to the reader that since [rik,-ff] ~ 0, the ground 
state of the system should be an eigenstate of nt, and therefore A^(k, k ) = for 
k 7^ k . This quantity then signifies the importance of terms beyond exchange 
energy. We have have tried an number of ansatzs for this object. In particular, if 
we ignore the dependence on the angle k.k , we may write a reasonable-looking 
formula for this object (k 7^ k ). 

7V(k,k') = ((nk) -nMk))((rv) -n^k')) (63) 

This has the attractive feature of conserving the number of particles and also 
the net momentum of the electrons ( -^(k, k ) sa 0, k A^(k, k ) w ) and 
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further since < rtk < 1 we expect — 1 < Af(k, k ) < 1 which is also obeyed. It 
has been pointed out that this object A^(k, k ) is singular in the BCS case 
. In particular, it is identically zero unless k = ±k . This means that for the 
BCS case which is a special case we may write, 

iV(k,k') = (7V) A(k)(5k,k'-<5k,-k') (64) 

Here A(k) is a function of order unity and (N) is the number of particles. 
Ideally, we would like a formula for this object that in the normal Fermi liq- 
uid regime reduces to an unremarkable function of |k|, |k | and k.k and ex- 
hibits a phase transition to the singular appearance as one crosses over into 
the superconducting regime. In fact it may be seen that the BCS correla- 
tion functions supports a net momentum in the sense of a quantum fluctuation 
{P'^)bcs = 2(A^) X^k^^^lk) 7^ 0. It will be shown later that such a formula 
is indeed possible in principle. Bare short-range interactions are not of much 
interest in continuum systems so we leave it at that (however Eq. ( ^ ) and 
Eq. ( ^ ) are meant to be universally valid ). 

3.3 Long-Range Interactions 

For long-range interactions, if one considers an approach similar to the one used 
for deriving the traditional RPA dielectric function, but this time including 
possible fluctuations in the momentum distribution, we find a result that is 
different form the usual RPA-dielectric function but one that reduces to it in 
the weak coupling limit. This derivation is found in Appendix D. We merely 
quote the final answer here: 

^eff{(l, u) = eg-RPA{(l, uj) ~ (-^) — ^. (65) 



Here, 



E 
k,k' 



and, 



PM.^)= (66) 
7V(k + q/2, k' + q/2) - 7V(k - q/2, k' + q/2) - JV(k + q/2, k' - q/2) + 7V(k - q/2, k' - q/2) 

(^ - ek-q/2 + ek+q/2) (t^ - Ck' -q/2 + ^k' +q/2) 

(67) 

1 , «q (no(k + q/2))-(rio(k^q/2)) 

^ k ^ ^ ek+q/2 + ek-q/2 

?k = ek-Ey(^k-q) (69) 

q#0 



The momentum distribution in such a system with long-range interactions may 
be written down as follows. 

fik = (k) (k) -f ( 1 - n^^ (k) ) (k) (70) 
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i^l(k) 

i^2(k) = 



1 



1 



1 



5b (k) 



q.l 



, l + 5B(k) 

^ SA(k) 



q) + k.q/m - eq)2 

1 - ^k+q 



y 

^ + k.q/m 

q,z ^ ^ ^/ 



5?(-q) 



5?(-q) 



3<(q) 



. k 



"-k4 



(^j(q) 



kq\2 



(71) 
(72) 
(73) 
(74) 

(75) 



This quantity (7i(q) is a residue of complex integration and is given as the 
frequency derivative of the polarization. 



?»(q) 



V 

Va 



d_ 



(76) 



here e^^^ is the principal part of e^f / . moreover cui are the zeros of this dielectric 
function. 



(77) 



As it stands the above equation is ill-defined since the dielectric function is 
complex. Interpreting uJi to be the zero of the principal part results in capturing 
only the collective mode. The particle-hole mode which is completely lost in this 
approach is very important. In order not to lose this mode we have to interpret 
the zeros in a special manner. We shall take the point of view that all positive 
energies are allowed as zeros but each comes with a weight corresponding to the 
strength of the dynamical structure factor at that energy. Define the weight to 
be 

1 



Ty(q, Lu) — Im 



.ee//(q,w-«0+) 
Therefore, sum over modes is now interpreted as, 



Y-./ N_/o duj W{q,uj)f{q,uj) 

Jo w^(q' ^) 



(78) 



(79) 



q,l 



It is worth emphasizing that Eq. ( |65| ) is valid for long range interactions where 
it is wrong to ignore the exchange correlation term. The one-particle properties 
that were derived in Eq.( p3) are valid for short range interactions only. The 
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interaction strength had to be weak so that one could ignore the exchange- 
correlation term in the short-range case. This was also needed in the long- 
range case to ensure that one could ignore the quadratic terms in the sea- 
displacements. However, in arriving at Eq.( ^5|) use was made of the form in 
Eq.( |5^) rather than Eq.( ^). There really is no a-priori reason to prefer one 
over the other. It is possible that one has to use a combination of the two (in 
fact, this has been shown in a preprint fl^). Using the second form exclusively 
(Eq.( leads to a dielectric function that is totally different from the RPA- 
dielectric function (we have not been able to write down an explicit formula 
for this), thereby casting doubt on the sanctity of the traditional methods of 
deriving dielectric functions. This issue has probably not been addressed fully by 
the many-body community even though a well-known text-book mentions some 
of these facts (Mahan|7| : end of chapter 4). One may argue that the latter 
form, namely in Eq.( ^6|) at least for Coulomb repulsion (1/r) in three spatial 
dimensions is negligible at high density. This is because in such a case we have 
two independent length scales kp^ and as and if kpas » 1, we may express 
all states as linear combinations of low lying excited states of the noninteracting 
system and then |k| « |k | « /cj- and < |q| ~ << kp, this means that 
|k - k'|2 >> {ag^y « g2 This means that Eq.( |5^) is negligible compared 
to Eq.( ^5|) . This ambiguity in operator ordering that leads to very different 
looking hamiltonians as in Eq.( ^) and Eq.( ^6|) suggests that this problem 
will persist in all cases where electron-electron interactions are present, in other 
words, in nearly all of many-body physics. This suggests to the authors that one 
must look elsewhere to find a more natural home for these techniques. We have 
some reason to believe that quantum electrodynamics is such a place. There, 
electron-electron interactions come about indirectly via-coupling to gauge fields 
where the operator ordering ambiguity is absent. Furthermore, a phenomenon 
important to condensed matter physics, namely magnetism, being primarily 
due to the spin of the electron, is taken into account naturally in the relativistic 
theory, since spin is a consequence of the Dirac equation^. In the next section, 
we study charge-conserving electron-hole systems that is a precursor to this 
more ambitious program of reworking gauge theories in the sea-displacement 
language. To conclude this section, we point out another definition of the sea- 
displacement that gives us the illusion of understanding the meaning of the 
square root. This is similar to Eq.( |l^) for bosons. 

-4k(q) = ni.(k-q/2)(l-nF(k4-q/2)) U - ^k-q/2+qi/2(qi)^k-q/2-Hqi/2(qi) 4-q/2Ck+q/2 

V qi#o / 

(80) 

The reason why this formula is not too illuminating is because a power series 
expansion around the unit operator and an iterative solution, yields just as in 
the Bose case, a result that amounts to ignoring the presence of the square root. 
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That is, 



Ak(q) =? np{k - q/2)(l - npik + q/2))cT _q/2Ck+q/2 (81) 

This is not the definition that is intended||l^. Instead, one must also consider 
regimes in which the square root is close to zero rather than close to unity. 
The exact definition of ^k(q) is the 'self-consistent' solution of the nonlinear 
operator equation Eq.( pO| ). Readers with superior mathematical skills may 
finally be able to solve this operator equation and map out its domain and 
range in Fock space and explain how this object acts on elements of its domain. 

3.4 Comparison With Other Truly Exactly Solved Models 

One of the other criticisms that has been expressed against the sea-boson 
method is that it is unable to reproduce Luttinger Liquid behavior. The authors 
have convinced themselves that the Calogero-Suthcrland Model(CSM) is in fact 
a Luttinger liquid. A Luttinger liquid is characterized by a lack of discontinuity 
in the momentum distribution, and the divergence (negative infinity) and con- 
tinuity of the slope of momentum distribution across the Fermi surface. The 
correlation functions of this model has been derived by Ha, Lesage et.al. and 
Forrester [Q. However, the correlation functions of this model depend on the 
use of a singular gauge transformation that forces a particular statistics ( in our 
case, fermionic ) on the correlation functions. The authors who have computed 
the correlation functions have used the natural statistics for the CSM which is 
in general, fractional. This makes comparison with our approach difficult. The 
CSM for statistical parameter /3 = 3 (in the notation of Lesage et.al. Q) is the 
simplest model that describes interacting fermions. However this choice makes 
the potential energy of the same order as the kinetic energy making the system 
strongly interacting. In our approach, the sea-bosons were written down in the 
plane wave basis with the tacit assumption that this is a good choice. It is easy 
to see that this choice is likely to be good only for weakly nonideal systems. The 
fermionic correlation functions for the weakly coupled CSM has yet to be writ- 
ten down. According to our expectations, this should result in a sharp Fermi 
surface for sufficiently weak repulsion. The Hubbard model in 1D[^ is another 
system that we have to compare our results with. It is exactly solvable only 
for half-filling. The correlation functions are hard to deduce in a closed form, 
but it is believed that this is a Luttinger liquid for arbitrarily weak repulsion. 
This latter result is the one we have the greatest difficulty in proving. In any 
event, in order for the community to accept the sea-boson technique,we have to 
somehow reproduce Luttinger liquid features. Work is in progress to do exactly 
this. It will be reported in a future publication. 
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4 Electron-Hole Systems : Exciton-Exciton In- 
teractions 



In this section, we study the exciton Green function using the sea-displacement 
technique. In fact, the biexciton Green function would be more interesting 
since it determines nonlinear optical response. This technique is ideally suited 
to study this quantity since exciton-exciton interactions are crucial in deter- 
mining the true nature of the bi-exciton. We find that exciton-exciton interac- 
tions are not of the two-body type as it is often assumed in the literature but 
rather more complicated. These interactions are mediated by other bosons that 
are present in the system, namely the intra-band particle hole excitations and 
the zero-momentum inter-band particle-hole excitation. All this will be made 
clear in the present section. Let us now move on to the details. As we have 
said before, a generalization is possible of these techniques to two-component 
charge-conserving electron-hole systems. One may define sea-displacement op- 
erators Akcr(qcr ) that include discrete internal degrees of freedom that may be 
metaphorically called spin. These objects are defined below. Since the definition 
in Eq.( [ij) was valid for both q = and q 7^ 0, we have, 

Aka(qa ) =nj^(k-q/2a)(l-nf(k + q/2a)) ^ 4-q/2agk+q/2.' (82) 

-y/^k— q/2(T 

We now make the following identifications. Let Ck be the operator that annihi- 
lates an electron from the conduction band and d^ y^ be the operator that creates 
a hole in valence band. Then we write, 

CkT = Ck (83) 

Cki = dlk (84) 

Then we have the following formulas for the electron-hole sea-displacement op- 
erators( again for both q = and q 7^ ) 

Aki(q T) = (l-<(-k-|-q/2))(l-7i|,(k+q/2)) ^ ^ c^-k+q/2Ck+q/2 

1 ~ ^-k+q/2'^-k+q/2 

(85) 

^kT(qi)=»F(k + q/2)n|.(k-q/2) ' =4-q/2^-k-q/2 (86) 

'Ck-q/2Ck-q/2 



AkT(q T) = "F(k - q/2)(l - n^(k + q/2)) =ci^_q^,Ck+q/2 (87) 

'^k-q/2'^k-q/2 



^ki(q i) = -n^(-k-q/2)(l-n^(-k+q/2))^=^^^^=dt_^_q/,d_k+q/2 

(88) 



i - ^lk+q/2^-k+q/2 
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If the system is undoped, we are tempted to set n|^(k) = n^(k) = 0. But 
this would be unwise. For such a choice would make all the sea-displacement 
operators but one, identically zero and this would lead to the conclusion that the 
exciton-exciton interactions are of the two-body type exclusively. We will soon 
argue that this underestimates the strength and importance of exciton-exciton 
interactions. The reason for this fallacy it seems is that we have to careful with 
the order in which we take limits. At any nonzero temperature, we expect that 
the k = state of the noninteracting system is always occupied even in a fully 
undoped system. Thermal fluctuations lead to a non-empty band. This in turn 
facilitates exciton-exciton interaction, since we may now contemplate intra-band 
particle-hole excitations competing for prominence with excitons and scattering 
off them and so on. In order to make all this more concrete, let us proceed as 
follows. We know that for a noninteracting Fermi system, 

n^(k) = 9{kF - |k|) n^(k) = e{kF - |k|) (89) 

If we go to the limit of an undoped system, fci? — * we find the following result, 

n^(k) = n^(k) = ^k,o (90) 

That is, only the zero momentum state is occupied, the rest are unoccupied. 
This choice enables us to have a scheme by which even in an undoped system, 
exciton-exciton interactions may be present and may contribute to the lincshape 
of the exciton. With this simplification we may write ( for both q = and q 7^ 



^ki(qT) = (l-'^k,q/2)(l-^k,-q/2) , =C?-k-Fq/2Ck-Fq/2 (91) 



k-Fq/2'^-k-Fq/2 



^kT(q i) = <5q,o<5k,o^=44 (92) 

^kT(q T) = 5k,q/2(l - 5k,-q/2)^^4cq (93) 
^ki(q i) = -5k,-q/2(l - <5k,q/2) ^-^ =dld^ (94) 

Let us now ascertain the commutation rules obeyed by these objects. Following 
the prescription outlined in Appendix B, we may write down the RPA-like 
commutation rules by replacing the right side of the approximate ones by the 
leading order results. If k 7^ ±q/2 and k ^ ±q /2 and q 7^ then, 

[Aki(q T), Ak'^(q' T)] = [Auiid T),4,^(q' T)] = ^k,k'^q,q' (95) 
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also, 

[Aki(OT),M(ot)] = o [Aki(0T),4^(0T)] = <5k,k' (96) 

[AoT(Oi),4^(Oi)]=l (97) 

[^q/2T(q T),^q'/2T(q T)] = [Aq/2T (Q T),^!,72t(^' t)] = ^q,q' (98) 

[^qm(q i),^q'/2i(q' i)] = [A^/21{CI i),^^/2^(q' i)] = 5q,q' (99) 

Finally, ^kt(0 T) = ^kj.(0 i) = . All other commutators involving any two 
of these sea-displacement operators are zero. In deriving these commutation 
rules, use has been made of the following approximate formulas. On the right 
hand side of commutation rules wc arc obliged to set the off diagonal terms to 
be identically zero and cjco = dgdo « 1. cj^Ck = c^^Ck « 0. It is worthwhile to 
verify some of these commutation rules explicitly. For example, 

[^ot(0 i), AI^{0 i)] « [cldl, doco] = clco - dodl = cJcq + d^do - 1 2 - 1 = 1 

(100) 

The rest arc reasonably straightforward. Using these facts, wc may write the 
following correspondence for the number conserving Fermi bilinears (the total 
number of electrons and holes commutes with these objects. See Appendix B 
for some hints as to how to derive these formulas). Here q 7^ 0, and we have 
singled out q = as a special case. If k ^ ±q/2 

4+q/2'=k-q/2 = ^ ^(l/2)(k+q/2)T(k + q/2 T)^(l/2)(k-q/2)T (k - q/2 T) 

+ ^ E <+q/2-qx/2i(qit)^k-q,/2i(-q + qiT) (101) 

qi74k+q/2 

^-k+q/2'^-k-q/2 = A ^|i/2)(k-q/2)i ("^ + q/2 i)^(l/2)(k+q/2)i ("k - q/2 i) 

+ ^ E <-q/2+q,/2i(qiT)^k+q,/2i(-q + qiT) (102) 

qi7t-k+q/2 

4c_q = ^4^0 ^-q/2T (-q T) « ^-q/2T (-q T) (103) 

4^0 = </2T (q t) \/^o « (q T) (104) 

= (i - £qd-q) ^ A^/^^ (-q i) « A^/^i (-q i) (105) 
dodl = At_^/2^(q i)(i - dld^)i « ALq/2i(q (106) 
Also we have for the number operators ( k ^ ). 

cjco = i - A ^ A[^^^^ (qi T)^q,/2T (qi T) - A 4^ (0 i) Apt (0 i) (107) 

qi#0 
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-qi/2i(qi T) (108) 

4do = i - A ^ ^^_q,/2i(qi i)^-qi/2i(qi i) - A 4^(0 i)yloT(0 i) (109) 
dUd-u = A 4/2i(-k i)A/2i(-k i) + A ^ 4+q,/2i(qi T)^k+qi/2i(qi T) 

qi7^-k 

(110) 

Now for the charge-conserving Fermi biHncar(tho total charge operator com- 
mutes with these objects). If k ±q/2 then, ( with q ) 

<^-k-q/2Ck-q/2 = (1 " rfLk-q/2^-k-q/2) ^^k^ (-q T) « -4ki(-q T) (HI) 

if in addition we have k ^ 

d_kCk = (i - d\d.^)^AuiiO T) « Aki(0 T) (112) 

if k = 

doco = Al^{0 I) ^[4^0- (Oi) (113) 

rfoC_q = A Al^iO i)A_q/2T(-q T)+A 4q,/2a(qi i)^-q/2-qi/2,i(-q+qi.T) 

(114) 

rf_qCo = -A Aj^(0 i)^q/2i(-q i)-A ^ 4i/2,t(^1 T)^q/2+qi/2a(-q+qi> T) 

qiT^O 

(115) 

Let us now verify some of these correspondences. For example we know that in 
the Fermi language, 

[4co,4c_q] = 4c_q (116) 

In the sea-displacement language we have, 

[4co,4c_q] = [yl_q/2t(-q T),4i/2t(^1 T)^qi/2T(qi T)] = ^-q/2T(-qT) ~ < 

qiT^O 

(117) 

as required. Similarly we have, 

[clco,doCo] = -doco (118) 
In the sea-displacement language we have, 

[clco,doCo] = -[4t(0 i)^0T(0 i),4T(0 i)] = -4t(0 i) ~ '^oco (119) 

Next, we write down the hamiltonian of the electron-hole system interacting via 
repulsion and attractive interactions. We show how excitons emerge naturally 
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from the formalism. We are also able to pin down the precise nature of exciton- 
exciton interactions. The total hamiltonian may be split into several parts. 
First is the kinetic energy plus the part of the potential energy that leads to the 
exciton. 

^0 = + ^9)4ck + E(^)^-k'^-k + Eg 4co 

E E 4w2'^-k'-,/2^-k'+q/2Ck-,/2 (120) 

q^O k±q/2/0-k'±q/2/0 

^^ = E^ E E 4+q/24'-q/2'^k'+q/2Ck-q/2 

q#0 k5<ik';k±q/2^0k'±q/2^0 

+ EI7 E E '^-k-q/2'^U'+q/2'^-k'-q/2'^-k+q/2 

q#0 k7ik';k±q/27i0k'±q/2^0 

_V-Wq V- (e) (e) _V-i^ V- (h) (h) .,21) 

2y ^ "k+q/2"'k-q/2 2^2^ "■-k-q/2"'-k+q/2 l^-^^-J 

qT^O kibq/2740 q#0 k±q/27i0 

-^-f.O = E ■y E (4cq+C^_qCo)4+q/2Ck-q/2 + E ^ E ('^-q^0+4c^q)c^lk+q/2'^-k-q/2 
q#0 kibq/2740 qT^O k±q/27t0 

- E ^ E (4+q/24'^qCk-q/2 + 4+q/2^^-q^0Ck-q/2) 

qT^O k±q/27t0 

-Ey E (^S^-k'+q/2'^-k'-c/2Cc, + cV-k'+q/2'^-k'-q/2^0) (122) 
qT^O -k'±q/27t0 

The part Hj when written out in terms of sea-displacements have four bosons 
in them. Whereas the term Hj^ has three bosons. The kinetic energy operator 
has only two bosons. Thus we may systematically regard the kinetic energy as 
being more important than Hj^ which in turn is more important than Hj. In 
the sea-displacement language, we may write upto additive constants, 

^0 = E(^K/2T(k t)^k/2T(k T) + E(^)<m(-k iMk/2i(-k i) 

+ E(£+^«)4i(0T)Aki(0T)- E ^4i(0T)Ak'i(0T) 

+ c^k(q)4i(qT)Aki(qT)- E ^ E 4i(q mk'i(q t) (123) 

k.qT^O k#k' qT^O 
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here, 



-.W = (I + ^. + - i;) + (1 - - (124) 

This hamiltonian has a very appealing form. The term containing Ak|(0 t) 
has been singled out since it highlights the exciton. This is nothing but the 
hamiltonian of the free exciton with center-of mass momentum equal to zero. 
The other hamiltonian involving Ak|(q t) corresponds to an exciton with non- 
zero center-of mass motion. If for example, we write, 

^ki(OT)=l]^/(k)6/(0) (125) 
I 

where <^/(k) is the Fourier transform of excitonic (hydrogenic) wavefunctions, 
then the hamiltonian may be recast in the diagonal form, 

(1^ + ^ ^^ki(o T)A'^(o T) =Y,Emb\mm 

(126) 

where Ei (0) are the energy levels of the exciton. we shall assume that a similar 
transformation has been performed for the full hamiltonian including center-of- 
mass motion. 

Ho = Y,EM)h\{o,)hM) - EgAl^{Q i)AoT(0 i) 
+ E(£-)4/2T(k T)A/2T(k T) + E(^)4/2i(-k i)^k/2i(-k i) (127) 

k^O k^O ^'"'^ 

Of particular interest in the term — Sg^J^(0 i)Ao|(0 J.). It suggests that the 
ground state of the system has one of these bosons present. This is consistent 
with the earlier observation since we have k = occupied by one electron and 
one hole. Since {Aq^{Q [)Y = it is clear that we can have only one of these 
bosons present. Let |0) be the vacuum of the fermions. Thus the Ao|(0 [) 
creates a solitron. The vacuum of the fermions may be written as, 

Ck|0) = d_k|0) = (128) 

It may be seen that this vacuum is the vacuum of all the bosons as well except 
^ot(0 i)- The action of Ao|(0 [) on |0) produces an electron-hole pair at k = 0. 
Since we have argued earlier that this corresponds to the ground state of the 
free theory, we have the following fact. The ground state of Hq is given by, 

|G)=AoT(Oi)|0) (129) 
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It may be seen quite easily that 

clco\G) = dldo\G) = 1\G) (130) 

The object ylk/2T (k T) annihilates a conductron. Finally the operator A_i(;/2j. (k i 
) annihilates a valeron. Now we would like to see how this evolves under the 
presence of the exciton-exciton interaction terms. A few words regarding these 
arc in order. Wc will sec here perhaps for the first time, the true nature of 
exciton-exciton interactions. They do not interact via simple two-body inter- 
actions rather the interaction is mediated by other bosons such as Ak/2t(k T) 
which arc intra-band particlc-holc excitations. Further wc have also coupling 
to the object Aq|(0 |) which is an inter-band zero momentum electron-hole ex- 
citation. Thus exciton-exciton interactions are mediated by these other bosons 
and the interactions arc rather more complex than the simple two-body variety. 
Let us write down precisely what they are. 

q^O k±q/2ji0 

+ (■*^/2T(«T) + '^U/2T<-'T)) J2 ''k+<./2-<,l/2i<''l'TMl'-'.l/2,i(-''+''l-T) 

qyiO k±q/2^0 qi;ik + q/2 

+ J2 (-'*q/2i<-''i'-^-j/2i(''i')'*a/2)(lc-q/2)i(-''+''/^'i'-*<l/2)(k+q/2)i(-''-'/2.i) 

q^O k±q/2^0 

+ 17 i-^lm^-" " - ^-'^m^" ■"Lq/2 + qi/2i<''l' T>'^k+qi/2,i(-« + «1. T) 

iTtO qi^-fc+q/2 

X] ^ ''(l/2)(k+q/2)T + "-^O.T i)'*(l/2)(k-q/2-q),i (k + q/2, T) 

q^O k±q/2^0 

X] 17 X] ^(l/2)(lc-q/2-q),i + "^0,T<"' i)'*(l/2)(k+q/2)T + <./2> T) 

q^O k±q/2^0 

Z-^ V /L-d (l/2)(k'-q/2),i (l/2)(k +q/2+q),i 

q^O -k'±q/2^0 

/ ' V / ' (l/2){k +q/2+q),i "^T (l/2)(k -q/2)l 

qTiO -k'±q/27!0 

"^X^ X Z]''(l/2)(k+q/2)-qi/2,i"'+'"''+''l'"-^-'»l/ = >i"'l " /2) (k-q/2 -q) , i 0^+"/^. 

q^O k±q/2=i0qi7!0 

"^Xi^ H X ^(l/2)(k-q/2-q),i"'+''^''"^-qi/2.if''^ ^ ( 1/2) (k+q/2) -qi /2 , i ^+''1 ■ T ) 

q^O k±q/25iOqi^O 

+ ^>^— >^ >^ / <-'''+''/2+<»l' T)'*q, /2 T*"! (-k'+q/2,T) 

/-^ V Z-^ /-^ (l/2)(k'-q/2)+qi/2,4. ^1/^.1 (l/2)(k +q/2+q),i 

17^0 -k'iq/2^0117^0 
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+ A > — \ > At (-k +q/2,T)At , (qi , ( - k' +q/2+qi , T ) 

V (l/2)(k'+q/2+q),i ^'^l ' ' > qi/2,T^''l" (l/2)(k' -q/2)+qi /2, i ' ^'^ ^ > 

q^O _k'±q/25iO<J15^0 

(131) 

"' = ''Ei7E E E ■"lv2)(kW2)T = • "^(V2)(k-q/2)T(''-<>/^>T)At^^^^^^,_^^^^^(k'-q/2.T)A^^^^,^^,^^^^,^^ 

17^0 fcyik' ''±<l/27^0k'±q/27^0 

+^'E^E E E E "U/2-q,/2a'-i"^-<.i/=a(-<.+<'i.»-;,,,,(,,_^,,,/>''-<'/^.T)^v2)(.'w^^ 

q^O k^k' k±q/2^0 k' ±q/27^0 H ^'<;+q/2 

+ ''E57E E E ^(l/2)(..W2)T'''+''/^'T'^(V2)(.c-q/2)T(W2,T)At,_^^^_^^^^^^(q2T)A^,_^^^^^^(^^^^ 

k^k' k±q/2j^0k'±q/2^0 

+ ^'E^E E E E E ■^U/2-qi/2,i<''^T>^''-'.l/ = .i<-''+'-"<'_q/2-q,/2a"'="V-q,/2.i(''+''-^ 

k;ik' k±q/2;i0 -|-q/2j^0qi^l' + q/2q2^k' -q/2 

q?^0 k^ik' I'±q/25i0k'±q/25i0 

+^'E^E E E E 

q^O kj^k' ''±q/27^0k'±q/27^0<ll?^-k+q/2 

17^0 k^ik' ''±q/27^0k'±q/25i0q27^-k'-q/2 

+ ^'E5^E E E E E "U/2+q,/2,i<''^"^-+<.l/ = .i'-+''-^K' "V+q,/2,i<''+''^ 

k^ik' l'±q/27^0k'±q/2^0ll?^-''+'5/2q2^-fc'-q/2 

'^'E 17 E ''(l/2)(k+q/2)T T)^(l/2)(k+q/2)T(>'+'>/2 "At^^^j^^_^^^j^ (k-q/2 T) A(i/2)(k-q/2)T (""O/^ T) 

q^O k±q/2^0 

-^'ES^ E E ■*a/2)(kW2)T'''+''/'"-*(V2)(k+q/2)T(''+'>/^TMt_^^^_^^^^^(q2T)A^_,/2-q2m 

q^O k±q/2^0q27!k-q/2 

"^'E^ E E 4+q/2-qi/2i<''l "-^k+q/a-qi/Sitqi T)At^^^^^^_^^^j^(k-q/2.T)A(i/2)(k-q/2),T(''-'</2T) 

q5^0 k±q/25i0qi^k+q/2 

-"'Ei^ E E E ^k+<./2-qi/2i(''lT>^k+q/2-qi/2i«.lT)At_^^^_^^^^^(q2T)A^_^/2-q2/2i('.2T) 

q5^0 k±q/2^0qi^k+q/2q2?!k-q/2 

"^'E 17 E "'(l/2)(k+q/2)i(-''— i)^(l/2)(k+q/2)i(-''-'>/^ ^'^(1/2) (k-q/2)i (-''+''/^ i)^(l/2)(k-q/2) i (-k+c/2 i) 

q^O -k±q/2j^0 

-^'ES^ E E ^(l/2)(k+q/2)i'-''-'/= ^'^(l/2)(k+q/2)i i)^Lq/2+q2/2iC2T)Ak-q/2+q2/2i<<>2T) 

q^O -k±q/25i0q2^-k+q/2 

-""'El^ E E 4+q/2 + qi/2i("'l "A^+^/2+qi/2i(<ll T)At^^^^^^__^^^j^(-k+q/2,i)A(i/2)(k-q/2)i(-''+<./2,i) 

q^O -k±q/27^0qi^-k-q/2 
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12 12 Yl 4+<,/2+c,i/2i(''lT)A,,+ <,/,+,^/2i(<>lTMLc,/2+<,2/2if''2T)A,,_^/,+,^/2i('.2T) 

q^O k±q/2^0qi^-fc-<l/2q2^-k+q/2 

(132) 

It is clear that this hamiltonian is rather more complex than the simple two-body 

variety. In particular, the leading contribution which is first order in A has three 
operators signifying a non-conservation of oxciton number. In one such process, an 
exciton recombines with a solitron to produce a conductron. Only cxcitons couple 
with external fields and we are able to infer the existence of the other excitations is 
inferred indirectly by studying the binding energy of biexcitons that are due to these 
exotic many-body processes. From the above form of the interaction terms we see that 
the ground state of Ho does not evolve with time under the action of the perturbation. 
However, if we first create an exciton (or two excitons) and then evolve the state under 
the action of the full hamiltonian, then the state does indeed evolve and produces all 
the effects that we expect, such as a possible exciton-line width, biexciton (bound state 
of two-excitons) and so on. Let us now try to verify these expectations rigorously. 

4.1 The Exciton Green Function 

In this section, we compute the exciton Green function as resulting from exciton- 
exciton interactions and compare the lifetime of the oxciton arising from exciton- 
exciton interactions with the lifetime from usual radiative recombination processes. 
Consider the initial state which corresponds to an exciton in internal state / and 
whose center of mass is moving with momentum Q. Then we may write, 

|/,Q)=6t(Q)|G) (133) 

Let us now examine how this state evolves with time. 

t (/f+Affj.,0 + ^^«/)|j.,Q) ^ * I i _ i A / dti/fj,o(ti) + {-i)2A2 / atiHj^o(tl> / dt2 H/,0 ('2 ) - « / dt i A j (tj ) ) i,+ (Q) | G> 

\ J Jo Jo Jo 

(134) 

There are many final states possible. We would like to compute the Green function, 

POO 

Gi{Q;uj)=i / dte'"*(/,Q|e-'*"|7,Q> (135) 



The exciton lineshape is given by plotting, 

L(w;Q) = -^7m(G/(Q;u;-i 0+)) (136) 
I 

We may evaluate the matrix element as, 

</,Q|e~* * (^0 + ^^f/,0 + ^^J^/)|j,Q) = j,-*ej(Q)t ^ ^_.^2^2^-Uj(Q)t ^^^,.•'^1 ^2 

X((5j(ki + qi/2 - Q/2; Q)<pj2 (ki - Q/2; -qi + Q)fj^ (^2 - Q/2; -qi + Q) ,5 j (kj + qi /2 - Q/2; Q) 
-i ee(qi)ti iej(Q)ti -"J-2(-qi+Q)tl / ^ i ee (qj ) "/j <Q "''1 '*2 . " " / (Q) t2 



Jo 







+ (-i)2^2e-"/(*3)' ^^(^)^<^*(ki-qi/2+Q/2;Q)v,j2(ki+Q/2; -qi +Q)¥> (''2 + Q/2; -qi +Q)»> j (kj -qi /2+Q/2; Q) 
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Jo Jo 

+ (-i)^^^e-"j(*3)t y~^(^)(I^)y*((l/2)Q-Qi;Q)v>j^((l/2)(Q-Qi)-qi;Q-Qi)v>*^((l/2)(Q-Qi)-q2;Q-Qi)¥.j((l/2)Q-Qi 

X f dti e-' «'.<«l)*le"-r«3>*le-"^2<-'3l+C!)tl /" ^* (Qi )t2 e"J-2 )*2 «3)*2 

Jo t/ 

+ (-i)2^2e-"-r(<3)t y^(^)(I^)y*(Qi-Q/2;Q)v>j^((l/2)(Qi-Q)+qi;0-Qi)y*^((l/2)(Qi-Q)+q2;Q-Qi)»>j(Qi-Q/2;Q) 
X / dti e-' «<^«3l)tl,ie_f(Q)ti^-«72(-«l+«)*l / ' dt^ e* «« «3l )*2 e"^2 (^-^l '*2 e""/ (C!)*2 

Jo Jo 

(137) 

We may rewrite the above equation retaining only the most singular parts as, 
(/Q|e""'^°+^^^°+^^^^'|/Q) = e~' * ^''•^^-i Fi{Q) t e~' ' ^''•'^^ « e"' * ('^^ 

(138) 

where, 

43*(ki + qi/2 - Q/2; QJcpj^ (''I - Q/2; -qi + Q)'!'/, ('^2 - Q/2i -qi + Q)<(>j(k2 + qi/2 - Q/2; Q) 



se(qi) + «72(Q - qi) - Ej-(Q) 

2 VJj(ki - qi/2 + Q/2; Q)-pj^ (kj + Q/2; -qi + Q)^*^{k2 + Q/2; -qi + Q)^j-(k2 - qi/2 + Q/2; Q) 



=E(^)(^) 



«fc(qi) + e/2 ~ "1^ ~ 'J'*-*^ 

*>3:(Ql - Q/2;Q)»>j-2((i)(<3i - Q) + qi; Q - Qi)-p*^ ((|)(Qi - Q) +q2iQ - Qi)v>i(Qi - Q/2; Q) 

«e(Ql) + «72 (Q - Ql) - '/(Q) 

v];(-Qi + Q/2;Q)v>72"i"*3 - Qi) - qi;Q - Qi)*.*^ ((1)(Q - Qi) - q2;Q - Qi)'Pj(-Qi + Q/2; Q) 



(139) 

Here ee(k) = k^/2me and ea(k) = /2mh- The lineshape of the exciton may be 
written as, 

L(0 a;^ = V I^iMQ)) 

In this article, we merely point out the feasibility of these computations. In our next 
article we intend to explore the practical consequences more thoroughly, specifically, 
the biexciton Green functions and nonlinear optical susceptibilities. 
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5 Conclusions 



In this article, we have laid down the ground work for a very promising new approach 
for understanding excitons in semiconductors and charge-conserving systems such as 
electrons and positrons. We have shown how to compute the cxciton-lincshapc that (if 
broadened) is due to also to nonradiative many-body processes. Wc have identified the 
elementary excitations in the system and precisely pointed out the role each of them 
play in the theory. The elementary entities in the two-component Fermi system may 
be thought of in two equivalent ways. One may consider them to be electrons and holes 
interacting via two-body attractive and repulsive interactions or we may consider the 
system to be made of excitons, solitrons, valerons and conductrons all interacting with 
each other by somewhat complicated, but purely local hamiltonians. The formidable 
technical challenges have been overcome and now a clear path has been mapped out 
for a systematic exploration of nonperturbative phenomena in many-body physics. 



6 Appendix A 

In this section, we prove the claims made in the first section namely that Xo is canon- 
ically conjugate to No- Let O = {bo — 

xo = i±t}ri^6^-i±t}r^d^^ 

2^ n 2^ n ^ ' 

n=l n=l 

n=l n=l 

n 

[d",No] = Y^6"-'^[d,No]6"'-^ (143) 

m=l 

But, 

[6,No] = i + d (144) 

Therefore, 

[()",iV„] =n d"~'(i + d) (145) 
It is a relatively simple matter to sum this series and we may write. 



[Xo, iVo] = ^ V ^ n 6"-i(i + o) + i V 6t"-i(i + &) 

n=l n=l 



= ^{i + d)-^ + ^(1 + 6^)^^ = 1 i (146) 

2' l + O 2"- 'l + Ot 

Therefore we find much to our relief, 

[Xo,iVo]=ii (147) 

Now we move on to fermions. 
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7 Appendix B 



Here we would like to prove some facts about fermions that were claimed in the 
main text. Writing down the commutation rules obeyed by v4k(q) involves knowing 
the precise meaning of the square root in the denominator. Since we have not been 
successful in pinning down its true meaning, we shall have to take the point of view 
that the meaning is uniquely fixed by demanding that the following commutation rules 
be obeyed. This is not very satisfying but the authors have exhausted their meagre 
capabilities. Let us consider the various possibilities, 
(a) nF(k + q/2) = 0, nF(k - q/2) = 1. 



V"k-q/2^k(q) = ( 1 - X] ^k-q/2+qi/2(qi)^k-q/2+qi/2(qi) 1 Ak(q) = 4_q/2Ck+q/2 
V qiT^O / 

(148) 

We have to convince ourselves that the left hand side obeys the same commutation 
rules as the right hand side. 

(b) npik + q/2) = 1, nF(k - q/2) = 0. 



-4|L(-q)V'^k+q/2 = ^I(-q' 



fi-E^Uq/ 



2+qi /2 (qi )^k+q/2+qi /2 (qi ) 



(c) nF(k + q/2) = 0, nF(k - q/2) = 0. 



4+q/2Ck-q/2 - X ^k+q/2-qi/2 

(d) nF(k + q/2) = 1, nF(k - q/2) = 1. 



(qi)^k-qi/2(-q + qi] 



' '=k-q/2'=k+q/2 



(149) 
(150) 



4+q/2''k-<l/2 = - X] ^k-q/2+qi/2(qi)^k+qi/2(-q + qi) 



Similarly, we have 

(o') nF(k' + q72) = 0, nF(k' - q72) = 1. 

v/^^v-q'/2^k'(q ) = 



■k'-q'/2 k +q /2 



(151) 



(152) 



and so on for the other cases. We would now like to write down some statements that 
would be analogous to commutation rules. Let us first define the following statements( 
it goes without saying that q 7^ and q ). 
SI : k + q/2 = k' +q'/2 
k + q/2 = k' - q72 
k- q/2 = k' + q72 
k- q/2 = k' - q72 

551 : k = k and q = q 

552 : k = k and q = — q 



S2 

S3 
S4 
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Consider the object : 

V«k-q/2^k(q)^%'_q'/2^k'(q') (153) 

We would now like to ascertain the meaning of this object when say, SI is true but 
S2, S3 and S4 are false. This state of affairs in symbolic logic is written as iSl A -'S2 A 
-.53 A -.54. Let us define, 

ASl = S1A ^52 A -nS'i A AS2 = S2 A -^Sl A ^53 A ^54 

ASS = 53 A -.51 A -.52 A -.54 ^54 = 54 A -.51 A -.52 A -.53 

AAO = ^51 A -^52 A -i53 A -^54 
(1) If ASl is true then we have, 

V»^k-q/2 Ak(q) ^n^' _q' /2^k' (q ) = y^%'_q'/2^k' (q )V"k-q/2^k(q) = (154) 



V'^k-q/2^k(q)A+,(-q')^nk/_^q//2 = (-q')yWk'+q'/2V'^k-q/2^k(q) = 

(155) 



v^^i:^^k(q) I 4'+q'/2-qi/2(*ii)^k'-qi/2(-q' +qi) 

\qi7^q',0 J 

= Ck-q/2Ck+q/2C^q/2Ck'-q'/2 = V"k-q/2^k_q' /2 (q-q ) I ^ ^ E ^k+q/2-qi/2(qi)^k+q/2-qi/2(qi) 

(156) 



E ^I'+q'/2-qi/2(^l)^k'-qi/2(-q' +qi) V^li^^k(q) 

^^qi5^q',0 / 

Ck+q/2Ck'-q'/24-q/2Ck+q/2 = " v''^k-q/2^k_q' /2 (q-q ) ( E ^k+q/2-qi /2 (qi)^k+q/2-qi/2 (qi) 

(157) 
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Vn^A^{ci){ E <'-c'/2+c,/2(qi)^k'+c,/2(-q'+qi) 1 =0 (158) 



E ^I'-,'/2+q./2(qi)^k'+c.i/2(-q' + V^ii^^k(q) = (159) 

i^qi#q',0 / 

^k(-q)V"k+q/2A;^, (-q')^%'+q'/2 = ^I' (-q')v'"k'w/2^I' (-q')y^vw7^ = ° 



(160) 



4(-q)V^^i^| E <'w/2-q./2(qi)^k'-q,/2(-q' + qi)| =0 (161) 



E 4'+q'/2-q,/2(^l)^k'-q,/2(-q' + qi ) J 4(-q)V^^l^ = (162) 
Vqi7^q',0 



'^ic(-0)V"''+'i/2 E ''I'-<,'/2+ci/2^''^'V+<,i/2(-'''+''l' j = ( E ■^k+<./2+ci/2(''l)'^>'+<./2+Cl/2(0l)J ^"''-"'^''k-c' /2 ' 

(163) 



/ ^ k -q /2+qi/2 k +qi/2 



At (-q)ynfc+q/2 = - I 1 - E '*k+q/2+qi /2 ''^k+q/2+qi /2 (H ) J y/rt^-.,/2\ 

\ qi^q.o / 



^qiTiq ,0 

(164) 

E 4+q/2-q,/2(qi)^k-q,/2(-q + qi)) I E +q' /2-q./2 (^1^)^^' -q, /2 (-q + q^) ) = 

l5^q,0 / \q2y^q,0 / 



(165) 



E <'W/2-q./2(^2)^k'-q,/2(-q'+q2)) ( E ^l+q/2-qi/2(qi)^k-q,/2(-q+qi)) =0 

2yiq,0 / \qi7^q,0 / 

(166) 

E ^k-q/2+qi/2(qiMk+q,/2(-q + qi)J E -q' /2+q. /2 +q. /2 ("q + q^) = 

l?^q,0 / \q2#q',0 / 



(167) 
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E <'-<,'/2+<,./2(^^)^k'+<.2/2(-q'+q2) ( E 4-c/2+,,/2(qiMk+q,/2(-q + qi) 

(168) 

Similarly, the reader may write down the corresponding relations when (2) AS2 is true 
(3) ASS is true and (4) when AS4 is true. If AAO is true then we have, 



[^k(q),A^'(q')] = [^k(q),4,(q')]=0 



If SSI is true then we have. 



1 - E ^k-q/2+qi/2(qi)^k-q/2+qi/2(qi) ^k(q) 



11^0 



(169) 



= (170) 



and for example. 



E ^k+q/2-qi/2(qi)^k-qi/2(-q + qi) j ( ^ ^k-q2/2(-q + q2)^k+q/2-q2/2(q2) j 
*vqi#q,0 / \q2y^q,0 / 



^k+q/2-qi/2 

'vqiy^O 



(qi)^k+q/2-qi/2(qi) 1- T At 



k-q/2-q2/2 



(q2)A 



k-q/2-q2/ 



2(q2) 



(171) 

Similarly, the reader can fill in the rest of the rules once the main techniques for 
deducing these rules have been laid down as we have done here. We see here that 
the objects Ak(q) and ^^(q) obey exact, closed commutation rules. This enables 
us to treat any theory involving fermions obeying simple fermion commutation rules 
in terms of a theory involving sea-displacements obeying these rather complicated- 
looking commutation rules. However, in all cases of practical interest we shall adopt an 
approximation that in the one-component Fermi system is equivalent to the RPA or its 
generalizations. In the multicomponent case, we may write down the following formula 
for the Fermi bilinears. This has been used to derive the corresponding formulas for 
the charge conserving electron-hole systems in the main text. If q 7^ then, 

4+q/2aCk-q/2a' = V"k+q/2a ^ka (-qO"' ) + ^L' (q'^) ^%-q/2a' 

+ E ^k+q/2-qi/2ai(qi«^)^k-qi/2ai(-q+qiO-')- ^ A^ (qifTl) A+qi/2a(-q+qi'Tl) 

qiT^q.OCTi qi7^q,0(Ti 

(172) 

4^Cka = njr(k0-)i-|- ^ ^k-qi/2ai(qiO')^k-qi/2ai(qiff)- ^ ^k+qi/2a (^l^^l )^k+qi/2a (qi^^l) 

(173) 



qiT^o.CTi 



qi7^o,CTi 
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7.1 Approximate Commutation Rules : The Random 
Phase Approximation 

The commutation rules presented above together with prescriptions such as Eq.( 
and Eq.( form a complete and closed system of rules. These rules however are 
far from simple to use. It is desirable to make some approximations in a systematic 
manner that enables us to simplify these rules. One such natural approximation is the 
RPA of Bohm and Pines [^]. The way to do this is to make the following assertion. 
Defn : The RPA is obtained by retaining only the lowest order sea-displacement 
terms in Eq.( and setting rik = n_F(k)i in Eq.( This means that in the exact 
definition in Eq.( |l^) we have to replace the number operator by the unit operator. 
When this is done, we have the RPA result for j4k(q). 

Ak(q) = npik - q/2)(l - Ufik + ci/2))cl_^^.^c^+^^2 (174) 



The igno ring of higher order terms amounts to the following (from Eq.( 150) and 

4+q/2Ck-q/2 = (175) 

when nF{k + q/2) = and nF(k — q/2) = OR when nF(k + q/2) — 1 and 
nF(k — q/2) — 1. Let us now use these simplifications to obtain a set of clo sed 
commutation rules valid in the RPA sense for the object Ak(q). From Eq.( 174) we 
have, 

[^k(q), A^' (q')] = npili - q/2)(l - ^^(k + q/2))7iF(k' - q72)(l - nF(k' + q72)) 

X [4-q/2'=k+q/2, C]^, _q/ /2Ck'+q'/2] 

= nF{k~ q/2)(l ^np{k + q/2))nF(k' - q72)(l - nF(k' + q72)) 

^ ('^k-q/2'^k'+q'/2'^k+q/2,k'-q'/2 ^ V -q' /2 '^''+l/2'^k-q/2,k' +q' /2 ) = ^ (^'''6) 

Now the right hand side is identically zero as we see from this argument. When 
k + q/2 = k' - q72 we have the factor (1 - nF(k + q/2))nF(k' - q /2) = 0. When 
k - q/2 = k' + q72 we have npCk - q/2)(l - ^^(k' + q72)) = 0. Therefore, 

[A^{ci),A^,{q)UpA = (177) 

Now let us compute, 

[4k(q),4, (q')] = UFik - q/2)(l - (k + q/2))nf (k - q72)(l ^ (k + q72)) 

^ [cLq/2'=k+q/2 , ^,C^ 

= nF(k - q/2)(l - nF(k + q/2))nF(k' - q72)(l - npi^l + q72)) 

^ ('^k-q/2'^k'-q'/2'^k+q/2,k'+q'/2 ^ '^I' +q' /2'^k+q/2 5k-q/2,k' -q' 72) (^"78) 



Because of Eq.( |175| ), the right hand side of the above equation ( Eq.( 178) ) is 
identically zero unless k = k and q = q . Therefore we may write, 

[Ak(q), (q )] = n^(k - q/2)(l - (k + q/2))5^ ^ (nk-q/2 - nk+q/2) (179) 
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If we retain the number operator as it is in the right hand side of Eq.( 179) we are 
dealing with the generalized RPA. The generalized RPA pays attention to possible 
fluctuations in the momentum distribution around a nonideal mean [^. However in 
the simple case we are obliged to set rik = uf (k)l- When this is done we find the 
following simple answer. 



[Ak(q), ^1^, (q')]flPA = nF(k - q/2)(l - npik + q/2))5 



Finally, we would like to address a rather important question, namely is the RPA 
as described above a controUed approximation ? Indeed, what really is a 'controlled 
approximation' ? And finally, does it matter whether some approximation is controlled 
or not ? In order to answer these questions it is important to first define what is 
meant by a controlled approximation. A definition that seems reasonable is as follows 
: An expansion in powers of a dimensionless parameter small compared to unity 
that is obtained by combining the various dimensionful parameters of the theory, is a 
controlled approximation. By this definition it is clear that the RPA is not a controlled 
approximation. Rather than expanding in powers of a dimensionless parameter it 
seems that we are expanding in powers of a dimensionless operator namely Ak(q) • 
For this to be justifiable, we have to show that this object is in some sense small 
compared to unity. This means that the matrix elements of this object have to be 
small compared to unity [^. This is possible only if we restrict our Hilbert space to 
be one that contains low lying excited stated of the noninteracting system. But this 
is really begging the question. We would like to know beforehand, given a type of 
interaction and its strength whether or not such an assumption is justified. In fact, in 
our earlier article |Q| we made some rather unfortunate remarks that may be forgiven 
since it was the first article in the series and we were going to fix the technical aspects 
later anyway. The formula for the 'sea-boson' given there in terms of the Fermi fields 
and justly criticized Cune and ApostolD is totally wrong despite its appealing form. 
The other remarks in the appendix|lj justifying the controlled nature of the RPA- 
approximation provided one restricts the functional form of the potential to obey the 
constraints outlined earlier [Q] so far seem alright. In any event, the moral of this 
discussion is that the RPA is a strange kind of approximation. It is not possible to 
make useful statements about when such an approximation breaks down. To say that 
RPA is valid if all states of the interacting system are expressible as linear combination 
of low lying states is as illuminating as saying : the best strategy to win a game of 
chess is to play it so that the opponent's king is placed under check and cannot move 
and the opponent cannot block or eliminate the check. 



8 Appendix C 

Consider the hamiltonian 

H = J2 + Y.W^ 4+q/2<'_q/2V+q/2Ck-q/2 (l8l) 
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here the exchange energy (minus the term quartic in the sea-displacements) has been 
absorbed into the kinetic energy. 

£k = ek-^^(nk_q) (182) 

Let nq(k) — c]^+q/2'^k-q/2 and no(k) = cJ^Ck Let us now write down the equation of 
motion for nq(k). Introduce a source, 

He.t{t) = 5^(?7e.t(qt) + U:,t{-qt)) ^ nq(k) (183) 

q k 

where Ue^tiqt) = e-'"*?7o(q). 

*^"q(k) = (ek-q/2-€k+q/2)<(k)+^ [Wq(k) , p^' ]p'_q/ + ^ 2^ Pq' KC^), P_q' ] 

qVo qVo 

+ iUe.Mt) + C/e*,«(-q't))[nq(k), (184) 

qVo 

["q(k),Pq'] = X]['^k+q/2Ck-q/2,C+,_^^,/2'=k'-q'/2] 
k' 

= X/'^k+<l/2'^'-q'/2'^k-q/2,k'+q'/2 ~ ''k' +q' /2'^''"'>/2^k+q/2,k' -q' /2 

k' k' 

« [no(k + q/2) - no(k - q/2)] (185) 

*^"q(k) = (ek-q/2 - ?k+q/2)nq(k) + y-(no(k + q/2) - n()(k - q/2))pq 

+ (t/e.t(-qt) + (7:,i(qt))(no(k + q/2) - no(k - q/2)) (186) 

Let us make a first pass at the computation of the dielectric function. Here, we make 
use of mean-field theory, that is, replace (no(k )pq) = (no(k ))(pq) 

«(k)) = (£k_q/2 - e\+q/2)«(k)) + ^((no(k + q/2)) - (no(k - q/2)))(p^) 

+ ((7e.t(-qt) + f/e*,t(qt))((no(k + q/2)) - (no(k- q/2))) (187) 
{n (k)) = ^ (no(k + q/2)) - (no(k - q/2)) ^ 

V U> — Ck-q/2 + fk+q/2 

+ (f/e.(-qt) + ^^.(qt)) <"°('- + -^/^)^-i"!('--^/^)^ (188) 

— ek-q/2 + ek+q/2 



This means. 



(P-q)=f/ext(qt)^^^ (189) 

Po(q,c.) = y (190) 
^ OJ — ek+q/2 + ek-q/2 

e(q,a;) = l-^Po(q,a;) (191) 
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Here, 



^From this and the fact that 

e,_HPA(q, a;) = = £(q, a;) (192) 

Next we would like to include fluctuations. Let us do this differently this time via the 
use of the BBGKY hierarchy. 

i^«(k)> = (?k-q/2 - €k+q/2)«(k)> + ^((no(k + q/2)) - (no(k - q/2)»(p^) 

+ ^(F2A(k + q/2, q) - F2A{k - q/2, q)) 
+ (f/e.t(-qt) + (7:,t(qi))(no(k + q/2) - no(k - q/2)) (193) 

F2A(k',q;t) = (no(k')p^) - (no(k')){p^) (194) 
F2(k';k,q;t) = (no(k')<(k)) - (no(k'))«(k)) (195) 
F2A(k',q;t) = ^F2(k';k,q;t) (196) 

k 

i^F2(k' ; k, q; t) = (?k-q/2-?k+q/2)i^2(k' ; k, q; t) + ^(N{k , k+q/2)-iV(k' , k-q/2))(p^) 

+ ^((no(k+q/2))-(no(k-q/2)))F2A(k',q;t)+([/e,«(-qf)+[/:,,(qt))(iV(k',k+q/2)-iV(k',k-q/2)) 

(197) 

Let us write, 

F2(k';k,q;t) = ?7e:.t(-qt)F2,a(k'; k, q) + ?7e*,t(qt)F2,i,(k'; k, q) (198) 

(pI) = [/e.t(-qi)(pq) + U:,Mt){p'^) (199) 

Also deflne, 

JV(k,k') = (no(k)no(k')) - (no(k))(no(k')) (200) 
u; F2,a(k';k,q) = (€k_q/2-ek+q/2)^^2,a(k'; k, q) + ^(iV(k' , k+q/2)-iV(k', k-q/2))(p;) 

+ ^((no(k+q/2))-{no(k-q/2)))F2"A(k',q) + (iV(k',k+q/2)-iV(k',k-q/2)) (201) 

-uj F2,6(k';k,q) = (ek_q/2-ek+q/2)J'2.b(k';k,q)+^(iV(k',k+q/2)-iV(k',k-q/2))(p;;) 

+ ^((no(k+q/2))-{no(k-q/2)))F|A(k',q) + (iV(k',k+q/2)-Ar(k',k-q/2)) (202) 

.(n u^) F'^ (U - V ^(k'.k + q/2)-iV(k,k-q/2) > 
e(q,w) i'2A(k ,<i)- TF / , 7, — ~ Tl 

V ^ Ul - ek-q/2 + ek+q/2 

^ y iV(k',k + q/2)-iV(k',k-q/2) ^^03) 

^ U — ek-q/2 + ek+q/2 

a; Cq(k) = (ek_,/2-ek+q/2)Cq(k) + ^((no(k+q/2))-(no(k-q/2)))(p;) + ^(F2"A(k+q/2)-F2"4(k-q/2)) 
+ ((no(k + q/2))-(no(k-q/2))) (204) 
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(no(k + q/2)) - (no(k - q/2)) , , 

V U) — ek-q/2 + ek+q/2 

(no(k + q/2))-{n„(k-q/2)) 



— ek-q/2 + ek+q/2 

«q (k + q/2, q) - F^^ (k - q/2, q) ^^05) 



V U) — ek-q/2 + ek+q/2 

After all this, it may be shown that the overall dielectric function including possible 
fluctuations in the momentum distribution is given by, 

/ \ / \ ^2(q,w) 

ee//(q, uj) = eg-HPA(q, w - T7 . (206) 

V eg-HPA (q,a;) 

Here, 

7V(k + q/2, k' + q/2) - iV(k - q/2, k' + q/2) - iV(k + q/2, k' - q/2) + iV(k - q/2, k - q/2) 

k,k' 



-e^ iV;.K + q/^,K +q/:^j - jV(K-q/Z,K + q/:^J - 7V + q/Z, K - cg-Z] 

ri I q, u>) = y 7 z z TT z : r 

^ (a; - ek_q/2 + ek+q/2 ) (t^ - e^' _q/2 + ^k' +q/2 ) 

'-,k 

(207) 

. Co , A - 1 -u (^o(k + q/2)) - (no(k - q/2)) 

ig-RPA(q,,iO) - i+ — y ; — (2U8) 
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